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Introduction 


The  metabolic  requirements  of  cancer  cells  are  different  than  those  of  normal  tissues  [1,2]. 
For  example,  the  liver  kinase  B1  (LKBl)  tumor  suppressor  protein  is  dysfunctional  in 
approximately  30%  of  non-small  cell  lung  carcinomas  (NSCLCs)  and  plays  a  critical  role  in  the 
cellular  response  to  energy  stress  [3].  Defects  in  such  metabolic  regulation  may  serve  as 
potential  points  for  therapeutic  intervention.  To  exploit  this  phenomenon  we  must  identify  the 
most  highly  activated  metabolic  fluxes  in  cancer  cells  as  a  function  of  their  genetic  background 
(e.g.  functional  status  of  the  LKBl -deficient  tumor  suppressor)  and  stress  level.  Systems-based 
approaches  are  required  to  quantitatively  describe  these  interconnected  pathways  in  sufficient 
detail  to  determine  optimal  enzyme  targets  under  conditions  of  stress  [4].  We  hypothesize  that 
combining  these  compounds  with  inhibitors  of  glutamine  metabolism  will  be  an  effective  therapy 
for  LKBl -deficient  lung  tumors.  The  objective  of  this  research  is  to  identify  metabolic 
susceptibilities  in  LKBl -deficient  NSCLC  using  systems-level  metabolic  flux  analysis. 
Subsequent  studies  in  Year  2  of  this  grant  will  aim  to  design  combinatorial  treatments  that 
simultaneously  induce  these  dependencies  and  target  compensatory  pathways. 

Body 

A  description  of  specific  results  associated  with  each  Task  scheduled  for  Year  1  is  included 
below.  Tasks  with  work  completed  or  in  progress  are  in  black  and  those  tasks  scheduled  for 
completion  in  Year  2  are  in  gray  font. 

Task  1.  Perform  systems-level  metabolic  analysis  of  isogenic  NSCLC  cell  lines  expressing  or 
lacking  functional  LKBl  using  isotope  tracers.  Lines  to  be  studied  include  vector  control 
(LKBl -deficient),  wild-type  (WT)  LKBl,  and  kinase  dead  (KD)  LKBl -expressing  A549 
cells.  Work  to  be  performed  by  Hui  Zhang  and  graduate  student  to  be  named  (Seth  Parker) 
(months  1-6) 

la.  Generate  C-labeling  data  in  the  above  cell  lines  using  glucose  and  glutamine  tracers 
under  conditions  of  nutrient  withdrawal  and  hypoxia  (months  1-4). 

lb.  Measure  nutrient  uptake  and  byproduct  secretion  rates  under  the  above  conditions 
(months  1-4). 

lc.  Validate  signaling  downstream  of  LKB  under  the  above  conditions  (months  1-4). 

ld.  Perform  computational  analyses  of  these  data  to  generate  comprehensive  flux  maps  of 
NSCLC  cells  +/-  LKBl  (months  3-6). 

We  have  completed  extensive  studies  comparing  A549  cells  engineered  with  an  empty 
pBABE  control  vector  (vec.)  versus  those  expressing  function  (wild-type)  LKBl  protein 
(+LKB1).  These  results  are  presented  in  Figure  1  below.  Although  the  cells  grow  identically  in 
culture  (Figure  lA),  addition  of  functional  LKBl  to  deficient  cells  changes  intracellular 
metabolic  pathway  activity.  Signaling  analysis  via  Western  blot  indicated  that  +LKB1  cells 
expressed  high  levels  of  the  tumor  suppressor  protein  (Figure  IB).  Culture  of  cells  in  the 
presence  of  [13C]glucose  (blue  circles)  or  [13C]glutamme  (green  circles)  and  mass  spectrometry 
analysis  (Figure  1C,  circles  represent  labeled/filled  or  unlabeled/open).  Quantitation  of  labeling 
in  TCA  intermediates  provides  a  relative  indication  of  metabolic  flux  from  tracers  to  specific 
metabolites.  For  example,  citric  acid  labeling  from  [13C]glucose  is  significantly  increased  in 
+LKB1  cells,  suggesting  that  flux  through  the  pyruvate  dehydrogenase  complex  (PDHC; 
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reaction  arrow  colored  red)  was  increased  when  functional  LKBl  was  present  in  NSCLC  cells. 
Consistent  with  this  expected  change,  we  observed  decreased  phosphorylation  of  acetyl- 
coenzyme  A  carboxylase  (pACC)  increased  but  increased  PDHC  phosphorylation  in  LKBl- 
deficient  A549s  (Figure  IE).  This  finding  is  consistent  with  a  role  of  hypoxia  inducible  factor-1 
(HIE- la)  in  reprogramming  metabolism  in  cells  lacking  functional  EKBl  [5].  Computational 
estimation  of  PDHC  flux  and  associated  confidence  intervals  from  sensitivity  analyses  further 
confirmed  this  finding  (Eigure  IE).  Metabolic  tracing  with  [13C]glutamine  indicated  there  was 
decreased  fully  labeled  (M+5  or  M5)  aKG  labeling  in  +LKB1  cells  (Eigure  IG).  Despite  this 
decrease  in  relative  glutamine  anaplerosis  (entry  of  carbon  into  the  TCA  cycle),  we  noted  a 
significant  increase  in  aKG  dehydrogenase  flux  (Eigure  IH)  and  increased  ATP-linked 
respiration  in  +LKB1  cells.  These  results  suggest  that  aKG  in  EKBl -deficient  cells  is  derived 
directly  from  glutamine  via  glutaminase  activity,  whereas  cells  with  a  WT  LKBl  protein  obtain 
aKG  predominantly  from  oxidative  isocitrate  dehydrogenase  (IDH)  flux. 
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Figure  1.  Metabolic  regulation  by  the  LKBl  tumor  suppressor.  A,  B)  Growth  rate  (A)  and 
LKBl  expression  (B)  in  A549  cells  expressing  an  empty  vector  (vec)  or  functional  LKBl 
(+LKB1).  C)  Atom  transition  map  depicting  isotopic  labeling  downstream  in  the  TCA  cycle 
from  [13C]glucose  (blue)  and  [13C]glutamine  (green).  D)  Citrate  labeling  from  glucose.  E) 
ACC  and  PDHC  phosphorylation  in  vec.  and  +LKB1  cells.  F)  PDHC  flux  calculated  via  MFA 
(error  bars  indicates  95%  confidence  interval).  G)  aKG  labeling  from  glutamine.  H)  aKGDH 
flux  calculated  via  MFA  (error  bars  indicates  95%  confldence  interval).  I)  ATP-linked  or 
uncoupled  respiration  in  vec.  and  +FKB1  cells.  Error  bars  indicate  standard  deviation  unless 
otherwise  noted. 


Additional  MFA  experiments  were  conducted  in  vec.  and  +LKB1  cells  under  conditions  of 
metabolic  stress,  including  hypoxia  or  matrix  detachment.  We  have  previously  demonstrated 
that  such  conditions  drive  metabolic  reprogramming  through  specific  glucose  and  glutamine- 
dependent  pathways  [6,  7].  While  NSCLC  expressing  or  lacking  LKBl  grew  at  similar  rates 
under  hypoxia  (Figure  2A),  [13C]glucose  tracing  indicated  that  PDHC  flux  was  differentially 
affected  under  this  condition  (Figure  2B),  with  +LKB1  cells  maintaining  higher  flux  through  this 
enzyme.  Nutrient/byproduct  uptake  and  secretion  were  not  significantly  different  (Figure  2C), 
but  the  extent  of  citrate  labeling  from  [13C]glutamine  was  reduced  in  -l-LKBl  cells  (Figure  2D), 
suggesting  these  cells  are  less  reliant  on  glutaminase  activity  than  vec.  control  cells.  Cells 
cultured  on  polyhydroxymethacrylate  (pHEMA)-coated  plates  are  unable  to  attach  to  matrix  and 
undergo  metabolic  stress  associated  with  decreased  PDH  flux  [6].  Expression  of  functional 
LKBl  indicates  that  these  cells  are  better  able  to  maintain  PDH  flux  (Figure  2E)  and  less  reliant 
on  glutaminase  (Figure  2F)  in  response  to  such  a  stress.  Time  course  analysis  of  AMPK,  ACC, 
and  PDH  phosphorylation  are  consistent  with  these  changes,  as  +LKB1  cells  have  increased 
pAMPK  and  pACC  but  decreased  pPDH  (Figure  2G).  These  results  indicate  LKBl  is  critical  for 
reprogramming  cell  metabolism  in  response  to  microenvironmental  stress. 
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Figure  2,  LKBl  mediates  metabolie  stress 
responses.  A-D)  Cell  growth  (A),  eitrate  labeling 
from  [13C]glueose  (B),  extraeellular  fluxes  (C),  and 
eitrate  labeling  from  [13C]glutamine  in  vee.  and 
+LKB 1  eells  under  normoxia  or  1  %  hypoxia.  E,F) 
Relative  PDH  (E)  and  glutaminase  (F)  fluxes  in 
vee.  and  +EKB 1  eells  24  hours  after  detaehment  on 
pHEMA  eoated  plates.  G)  Phosphorylation  of 
AMPK,  ACC,  and  PDH  over  24  hours  in  vee.  (-)  or 
+EKB1  (+)  eells  over  24  hours  after  plaeement  on 
pHEMA  eoated  plates.  Error  bars  indieate  SD. 
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Task  2.  Observe  the  effects  of  biguanide  treatment  on  NSCLC  cell  lines  at  physiologically 
relevant  doses.  NSCLC  cell  lines  to  be  studied  include  the  above  A549  cells  and  a  similar 
panel  of  H460  cells.  Work  to  be  performed  by  Hui  Zhang  and  graduate  student  Seth  Parker 
(months  6-10) 

2a.  Perform  '^C-tracing  experiments  above  upon  treating  cells  with  decreasing  doses  of 
metformin  and  phenformin  (months  6-10) 

2b.  Computationally  resolve  intracellular  fluxes  using  the  above  data  (months  6-10) 

2c.  Quantify  cell  survival  in  response  to  nutrient  deprivation  and  hypoxia  at  the  relevant  does 
of  biguanides  (months  6-10) 


We  have  conducted  MFA  studies  tracing  the  reliance  of  NSCLC  cells  on  glutamine  for  AcCoA 
generation  in  response  to  phenformin  treatment  and  in  the  context  of  cells  that  are  deficient 
(vec.)  or  proficient  (-I-LKB1)  in  the  LKBl  tumor  suppressor.  [13C]glutamine  labeling  patterns 
indicative  of  oxidative  glutaminolysis  (blue)  or  reductive  carboxylation  (red)  are  depicted  in 
Figure  3A.  Phenformin  drives  a  pronounced  increase  in  glutaminase  flux,  as  indicated  by  the 
extent  of  fully  labeled  (M5)  glutamate  from  [13C]glutamine  tracing  (Figure  3B).  Importantly, 
+LKB 1  cells  more  proficient  in  synthesizing  asparate  from  glutamine  in  the  presence  of 
phenformin,  as  evidenced  by  the  decreased  MO  and  increased  M3  labeling  in  -l-LKBl-l-phen  cells 
(Figure  3C).  Computational  MFA  was  conducted  on  palmitate  labeling  and  demonstrated  that 
+LKB1  cells  maintained  slightly  higher  glucose  oxidation  and  glutamine  conversion  to  lipids 
(Figure  3D),  providing  evidence  that  they  are  better  able  to  survive  and  proliferate  in  the 
presence  of  phenformin.  In  fact,  +LKBI  cells  were  able  to  proliferate  (albeit  at  a  slow  rate)  with 
50  pM  phenformin  in  the  culture  medium,  whereas  cells  lacking  functional  LKBl  did  not 
proliferate  at  all  (Figure  3D).  These  results  support  our  hypothesis  that  phenformin  treatment 
and  glutamine  metabolism  may  be  effective  co-targets  to  slow  LKBl -deficient  tumor  growth. 
Similar  results  were  also  observed  using  H460  NSCLC  cells  lacking  or  expressing  functional 
LKBl  (not  shown). 
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Figure  3.  Phenformin  reprograms  glutamine  metabolism.  A)  Atom  transition  map  depicting 
oxidative  (blue)  and  reductive  (red)  metabolic  labeling  from  [13C]glutamine.  B)  Glutamate 
labeling  from  [13C]glutamine.  C)  Aspartate  labeling  from  glutamine.  D)  MFA  modeling 
quantifying  the  contribution  of  glucose  and  glutamine  to  lipid  carbon  (error  bars  indicates  95% 
confidence  interval).  E)  Growth  rate  of  vec.  and  +LKB1  cells  in  the  absence  or  presence  of  50 
pM  phenformin.  Error  bars  indicate  standard  deviation  unless  otherwise  noted. 


Task  3.  Identify  the  speeifie  enzymes  which  catalyze  glutamine  anaplerosis  in  response  to  stress 

Work  to  be  performed  by  Hui  Zhang  and  graduate  student  Seth  Parker  (months  10-14). 

3a.  Quantitatively  compare  the  extent  of  glutamine  versus  glutamate  uptake  in  NSCLC  cells 
under  stress  (months  10-11) 

3b.  Quantify  and  compare  fluxes  through  glutaminase,  nucleotide  synthesis,  and  hexosamine 
synthesis  using  [y-^^N] glutamine  to  determine  the  enzyme(s)  which  convert  glutamine  to 
glutamate  (months  11-12) 

3c.  Quantify  and  compare  the  flux  of  glutamate  to  a-ketoglutarate  as  catalyzed  by 

aminotransferases  and/or  glutamate  dehydrogenase  using  [a-^^N]glutamine  (months  12- 
13) 

3d.  Functionally  validate  the  effects  of  metabolic  inhibitors  of  each  enzyme  using  MFA  and 
treatments  with  asparaginase,  968,  BPTES,  azaserine,  DON,  aminooxyacetate,  or  EGCG 
(months  12-14). 
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Specific  glutamine  and  glutamate  tracing  studies  using  [13C]  and  [15N]  isotopes  continue  to  be 
executed.  Given  the  results  in  Tasks  1  and  2  that  suggest  glutaminase  activity  is  elevated  in 
LKBl -deficient  NSCLC  cells  we  have  accelerated  experiments  that  utilize  inhibitors  of  this 
enzyme.  BPTES  and  968  are  published  inhibitors  of  glutaminase,  but  the  former  compound  is  a 
particularly  good  candidate,  as  recent  studies  (personal  communication  with  Chi  Dang,  UPenn) 
indicate  BPTES  is  much  more  effective  at  inhibiting  glutaminase  in  culture  and  in  vivo. 
Preliminary  titration  studies  quantifying  cell  growth  inhibition  by  BPTES  indicate  that  vec.  cells 
are  not  more  sensitive  to  BPTES  than  +LKB1  cells  (Eigure  4).  However,  the  sensitivity  may  be 
different  in  the  context  of  phenformin  treatment  or  within  in  vivo  microenvironments,  as 
preliminary  studies  that  are  being  repeated  suggest  that  the  recovery  of  vec.  cells  in  response  to 
phenformin  treatment/washout  is  compromised  relative  to  +EKB1  cells  and  the  glutaminase- 
dependence  is  exacerbated  under  such  conditions. 


Figure  4,  Cell  growth  inhibition  by  BPTES.  Cells 
were  growth  in  varying  eoneentrations  of  the 
glutaminase  inhibitor  BPTES  and  quantified  via 
plate-eounting.  While  BPTES  redueed  the  growth 
of  A549  eells,  under  non-stress  eonditions  no 
seleetive  sensitivity  was  observed.  Error  bars 
indieate  standard  deviation. 


Task  4.  Evaluate  the  efficacy  of  combinatorial  treatments  in  EKBl -proficient  and  deficient 
NSCLC  cells.  Work  to  be  performed  by  Hui  Zhang  and  graduate  student  TBN  (months  14- 
16). 

4a.  Perform  cell  survival  studies  by  treating  cells  with  combinations  of  biguanides  and 
inhibitors  of  glutamine  metabolism  (months  14-16) 

4b.  Analyze  the  metabolic  phenotype  of  cells  upon  treatment  with  drug  combinations 
(months  14-16) 

4c.  Identify  the  optimal  target/drug  to  be  used  in  combination  therapy  (months  15-16) 

Task  5.  In  vivo  testing  of  biguanide  and  glutamine  metabolism  inhibitors  in  xenograft  models  of 
LKBl -proficient  and  deficient  NSCLC.  Work  to  be  performed  by  Rob  Svensson  at  Salk  and 
Hui  Zhang  and  graduate  student  TBN 

5a.  Obtain  lACUC  and  ACURO  approval  for  combination  therapies  in  xenograft  mice 
(months  12-15) 

lACUC  and  ACURO  approval  have  been  granted  for  in  vivo  xenograft  studies,  which  will 
commence  in  the  next  3  months. 

5b.  Initiate  xenograft  experiments  with  A549  cell  line  panel  (control,  WT  LKBl,  KD  LKBl) 
and  treat  with  vehicle.  Gin  inhibitor,  phenformin,  or  phenformin  +  Gin  inhibitor  (months 
15-19) 

5c.  Conduct  xenograft/combo  therapy  experiments  with  H460  cell  line  panel  (months  17-21) 
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5d.  Evaluate  the  tumor  burden  and  eell  signaling  profile  of  treated/untreated  tumors  through 
lueiferase  and  histology  analyses  (months  18-23) 

5e.  Quantify  metabolite  profiles  in  plasma,  tumors,  and  adjaeent  normal  tissue  using  GC/MS 
analysis  (months  18-24) 

This  grant  partially  supported  work  on  related  13C  MFA  projeets  aimed  at  studying  the 
metabolic  response  of  cells  to  hypoxia  in  the  presence  or  absence  of  mutant  IDHl  and  IDH2. 
Additional  studies  developed  a  new  approach  for  quantifying  metabolic  sources  of  NADPH 
regeneration  in  the  mitochondria  and  cytosol.  Both  of  which  resulted  in  publications  [8,  9]. 

Key  Research  Accomplishments 

•  Identified  key  differences  in  metabolic  activity  in  NSCLC  cells  lacking  functional  LKBl 

•  Demonstrated  key  stress  conditions  that  exacerbate  the  defect  caused  by  LKBl  loss 

•  Uncovered  a  mechanistic  link  between  LKBl  activity  and  maintenance  of  PDH  flux 

•  Characterized  metabolic  reprogramming  in  response  to  phenformin 

•  Explored  sensitivity  of  LKBl -proficient  and  deficient  cells  to  glutaminase  inhibition 
under  nutrient  replete  conditions 

•  Demonstrated  the  increased  dependence  of  hypoxic  cells  on  both  oxidative  and  reductive 
glutamine  metabolism 

•  Developed  a  tool  for  quantifying  compartment-specific  NADPH  metabolism  in  the 
cytosol  and  mitochondria 

Reportable  outcomes 

Research  presentations: 

9/2013  FSU  Biochemistry  Seminar,  Tallahassee,  FA 

10/2013  Agios  Pharmaceuticals  MiniSymposium,  Cambridge,  MA 

3/2014  Keystone  Tumor  Metabolism  meeting,  Whistler,  Canada 

3/2014  Vesalius  Research  Center,  Leuven,  Belgium 

3/2014  Luxembourg  Centre  for  Systems  Biology 

4/20 1 4  AACR  National  Meeting 

4/2014  use  Biomedical  Engineering  Department  Seminar 

Two  research  manuscripts  were  published  with  partial  support  from  this  grant  [8,  9]. 

Conclusions 

Our  MFA  results  have  generated  a  “metabolic  roadmap”  highlighting  the  distinct  metabolic 
changes  that  are  caused  by  loss  of  the  LKBl  tumor  suppressor  in  NSCLC.  These  data  indicate 
that  LKBl -deficient  cells  have  a  propensity  for  using  glutamine  to  fuel  bioenergetic  metabolism 
in  the  TCA  cycle.  We  have  moved  beyond  this  basic  characterization  and  also  investigated  how 
metabolism  changes  in  LKBl -deficient  and  proficient  NSCLC  cells  when  they  experience 
stresses  that  affect  the  energy  state  of  cells,  including  hypoxia,  matrix  detachment,  and  treatment 
with  the  biguanide  phenformin.  These  studies  indicate  that  metabolic  stress  further  exacerbates 
the  metabolic  differences  between  NSCLC  cells  that  express  or  lack  functional  LKBl .  Of 
particular  interest  is  the  reliance  of  LKBl -deficient  cells  on  flux  through  the  glutaminase 
enzyme.  Inhibitors  are  under  investigation  in  the  clinic  in  some  cancers,  and  we  are  now 
exploring  the  efficacy  of  such  treatments  in  LKBl -deficient  cells.  While  marginal  differences  in 
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response  have  been  observed  thus  far,  these  experiments  have  relied  thus  far  on  eell  culture 
growth  only.  We  hypothesize  that  the  in  vivo  microenvironment  may  drive  LKBl -deficient  cells 
to  be  more  sensitive  to  metabolic  stress  or  combination  treatments.  Studies  using  xenograft 
models  of  lung  tumor  growth  will  commence  in  Year  2  of  this  grant. 

Ultimately,  these  results  can  make  a  rapid,  clinical  impact  by  outlining  specific  genetic 
phenotypes  in  NSCLC  that  may  benefit  from  combinatorial  treatments  with  biguandies  such  as 
metformin  (clinically  used  for  diabetes  already)  or  phenformin  and  glutaminase  inhibitors. 
Several  compounds  that  inhibit  GLS  are  currently  under  investigation  at  the  pre-clinical  and 
Phase  I  stages  (Calithera).  However,  metabolic  treatments  are  likely  to  most  effective  when  the 
specific  regulatory  or  enzymatic  defects  of  a  given  cancer  are  exploited  as  outlined  here.  Thus, 
our  combinatorial  therapeutic  approach  may  provide  a  more  efficacious  means  of  exploiting 
cancer  metabolism. 

References 

1.  DeBerardinis,  R.J.,  et  ah.  The  biology  of  cancer:  metabolic  reprogramming  fuels  cell 
growth  and  proliferation.  Cell  Metab,  2008.  7(1);  p.  11-20. 

2.  Vander  Heiden,  M.G.,  L.C.  Cantley,  and  C.B.  Thompson,  Understanding  the  Warburg 
effect:  the  metabolic  requirements  of  cell  proliferation.  Science,  2009.  324(5930):  p. 
1029-33. 

3.  Shackelford,  D.B.  and  R.J.  Shaw,  The  LKBl-AMPK pathway:  metabolism  and  growth 
control  in  tumour  suppression.  Nat  Rev  Cancer,  2009.  9(8);  p.  563-75. 

4.  Zamboni,  N.,  (13)C  metabolic  flux  analysis  in  complex  systems.  Curr  Opin  Biotechnol, 

2010. 

5.  Shackelford,  D.B.,  et  ah,  mTOR  and HIF-1  alpha-mediated  tumor  metabolism  in  an  LKBl 
mouse  model  of  Peutz-Jeghers  syndrome.  Proc  Natl  Acad  Sci  USA,  2009.  106(27);  p. 
11137-42. 

6.  Grassian,  A.R.,  et  ah,  Erk  regulation  of pyruvate  dehydrogenase  flux  through  PDK4 
modulates  cell  proliferation.  Genes  Dev,  2011.  25(16):  p.  1716-33. 

7.  Metallo,  C.M.,  et  ah.  Reductive  glutamine  metabolism  by  IDHl  mediates  lipogenesis 
under  hypoxia.  Nature,  2012.  481(7381);  p.  380-4. 

8.  Lewis,  C.A.,  et  ah.  Tracing  Compartmentalized NADPH Metabolism  in  the  Cytosol  and 
Mitochondria  of  Mammalian  Cells.  Mol  Cell,  2014. 

9.  Grassian,  A.R.,  et  ah,  IDHl  Mutations  Alter  Citric  Acid  Cycle  Metabolism  and  Increase 
Dependence  on  Oxidative  Mitochondrial  Metabolism.  Cancer  Res,  2014.  74(12):  p.  3317- 
31. 


11 


Please  cite  this  article  in  press  as;  Lewis  et  ai.,  Tracing  Compartmentaiized  NADPH  Metabolism  in  the  Cytosol  and  Mitochondria  of  Mammalian  Cells, 
Molecular  Cell  (2014),  http://dx.doi.Org/10.1016/j.molcel.2014.05.008 

Molecular  Cell 

Article 


CelPress 


Tracing  Compartmentalized  NADPH  Metabolism 
in  the  Cytosol  and  Mitochondria  of  Mammalian  Cells 

Caroline  A.  Lewis, Seth  J.  Parker, 2-®  Brian  P.  Fiske,^  Douglas  McCloskey,^  Dan  Y.  Gui,’  Courtney  R.  Green,^ 
Natalie  I.  Vokes,  Adam  M.  Feist, ^  Matthew  G.  Vender  Heiden,^’®*  and  Christian  M.  Metallo^’"^* 

■'The  Koch  Institute  for  Integrative  Cancer  Research  at  Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139,  USA 
^Department  of  Bioengineering,  University  of  California,  San  Diego,  La  Jolla,  CA  92093,  USA 
'^Dana-Farber  Cancer  Institute,  Boston,  MA  02115,  USA 

''Institute  of  Engineering  and  Medicine,  University  of  California,  San  Diego,  La  Jolla,  CA  92093,  USA 
®These  authors  contributed  equally  to  this  work 

‘Correspondence:  mvh@mit.edu  (M.G.V.H.),  cmetallo@ucsd.edu  (C.M.M.) 
http://dx.doi.Org/10.1016/j.molcel.2014.05.008 


SUMMARY 

Eukaryotic  cells  compartmentalize  biochemical  pro¬ 
cesses  in  different  organelles,  often  relying  on  meta¬ 
bolic  cycles  to  shuttle  reducing  equivalents  across 
intracellular  membranes.  NADPH  serves  as  the  elec¬ 
tron  carrier  for  the  maintenance  of  redox  homeosta¬ 
sis  and  reductive  biosynthesis,  with  separate  cyto¬ 
solic  and  mitochondrial  pools  providing  reducing 
power  in  each  respective  location.  This  cellular  orga¬ 
nization  is  critical  for  numerous  functions  but  compli¬ 
cates  analysis  of  metabolic  pathways  using  available 
methods.  Here  we  develop  an  approach  to  resolve 
NADP(H)-dependent  pathways  present  within  both 
the  cytosol  and  the  mitochondria.  By  tracing  hy¬ 
drogen  in  compartmentalized  reactions  that  use 
NADPH  as  a  cofactor,  including  the  production  of 
2-hydroxyglutarate  by  mutant  isocitrate  dehydroge¬ 
nase  enzymes,  we  can  observe  metabolic  pathway 
activity  in  these  distinct  cellular  compartments. 
Using  this  system  we  determine  the  direction  of 
serine/glycine  interconversion  within  the  mitochon¬ 
dria  and  cytosol,  highlighting  the  ability  of  this 
approach  to  resolve  compartmentalized  reactions 
in  intact  cells. 


INTRODUCTION 

One  of  the  defining  characteristics  of  eukaryotic  celi  metaboiism 
is  the  compartmentaiization  of  reactions  in  different  organelies. 
Although  coordination  of  metabolic  flux  across  organelles  is  crit¬ 
ical  for  cell  physiology,  the  inability  to  distinctly  observe  identical 
reactions  present  in  more  than  one  subcellular  location  has  been 
a  major  barrier  to  understanding  cell  metabolism.  Many  of  these 
compartmentalized  reactions  are  oxidation/reduction  (redox)  re¬ 
actions  that  utilize  pyridine  nucleotide-based  cofactors  to  trans¬ 
fer  electrons  between  metabolites  to  support  biosynthesis, 
redox  homeostasis,  signal  transduction,  and  ATP  generation 
(Poliak  et  al.,  2007a).  For  instance,  reduction  of  NAD'‘  to  NADH 


captures  energy  from  catabolic  reactions  to  drive  ATP  synthesis 
through  mitochondrial  oxidative  phosphorylation,  while  NADPH 
is  regenerated  via  a  different  set  of  reactions  to  maintain  reduced 
glutathione  (GSH)  pools  and  support  reductive  biosynthesis 
(Lunt  and  Vander  Heiden,  2011).  As  such,  NADPH  has  been  hy¬ 
pothesized  to  be  limiting  for  proliferation,  lipid  biosynthesis,  and 
survival  in  response  to  cell  stress  (Diehn  et  al.,  2009;  Jeon  et  al., 
2012;  Jiang  et  al.,  2013;  Schafer  et  al.,  2009).  These  compart¬ 
mentalized  metabolic  processes  impact  numerous  cell  and 
tissue  functions;  therefore,  understanding  how  biochemical  net¬ 
works  function  across  compartments  is  necessary  to  determine 
how  metabolism  contributes  to  disease  pathologies. 

The  pool  of  NADP(H)  in  cells  is  small  relative  to  flux  through 
pathways  that  utilize  this  cofactor  (Poliak  et  al.,  2007a).  Thus, 
interconversion  between  the  oxidized  and  reduced  states  must 
be  coupled  across  all  reactions  involving  this  cofactor,  and 
changes  in  abundance  may  not  be  informative  for  assessing 
the  use  of  NADPH  in  a  particular  pathway.  Neither  NAD(H)  nor 
NADP(H)  is  known  to  be  transported  across  intracellular  mem¬ 
branes  (Nikiforov  et  al.,  201 1 ;  Poliak  et  al.,  2007b),  and  multistep 
shuttles  involving  compartmentalized  redox  reactions  are  used 
to  transfer  electrons  between  the  mitochondria  and  cytosol  (Bis- 
sell  et  al.,  1976;  LaNoue  et  al.,  1974;  LaNoue  and  Schoolwerth, 
1979).  This  organization  facilitates  the  maintenance  of  different 
NADPH/NADP’^  ratios  in  each  subcellular  location  and  allows 
for  the  execution  of  compartment-specific  metabolic  processes. 
Classically,  cytosolic  NADPH  is  thought  to  be  regenerated  pri¬ 
marily  via  the  oxidative  pentose  phosphate  pathway  (PPP) 
(Lunt  and  Vander  Heiden,  201 1 ;  Poliak  et  al.,  2007a).  Other  po¬ 
tential  sources  of  cytoplasmic  NADPH  exist  in  mammalian  cells, 
including  reactions  catalyzed  by  specific  isozymes  of  isocitrate 
dehydrogenase  (IDH),  malic  enzyme  (ME),  aldehyde  dehydroge¬ 
nase  (ALDH),  and  methylene  tetrahydrofolate  dehydrogenase 
(MTHFD)  (Poliak  et  al.,  2007a;  Tibbetts  and  Appling,  201 0).  How¬ 
ever,  isoforms  of  several  of  these  enzymes  also  catalyze  iden¬ 
tical  reactions  in  the  mitochondria  and  can  potentially  transfer 
reducing  equivalents  between  the  mitochondria  and  the  cytosol. 
For  example,  the  reductive  carboxylation  of  alpha-ketoglutarate 
(otKG)  to  isocitrate  by  IDH2  consumes  mitochondrial  NADPH, 
with  citrate/isocitrate  subsequently  transported  to  the  cytosol 
where  it  can  be  oxidized  by  IDH1  to  produce  cytosolic  NADPH 
(Sazanov  and  Jackson,  1994;  Wise  et  al.,  2011).  Theoretically, 
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Figure  1.  Use  of  Glucose  to  Label  Cytosolic  NADPH 

(A)  Atom-transition  map  depicting  a  model  of  deuterium  transfer  from  [3-^H]glucose  through  glycolysis  and  the  pentose  phosphate  pathway.  Open  large  circles 
represent  carbon,  and  small  red  circles  indicate  deuterium  label  from  [3-^H]glucose.  Where  measured,  enrichment  of  Ml  isotopomer  (%)  for  glycolytic  in¬ 
termediates  in  parental  H1299  cells  is  shown. 

(B)  Labeling  of  NADPH  from  [3-^H]glucose  in  parental  HI 299  cells  overtime. 

(C)  Saturated  fatty  acid  labeling  (myristate;  Cl  4:0,  palmitate;  Cl  6:0  and  stearate  C18:0)  from  [3-^H]glucose  in  parental  HI  299  cells  following  incubation  for  72  hr. 

(D)  Cholesterol  labeling  from  [3-^H]glucose  in  parental  H1299  cells  cultured  for  72  hr. 

(E)  Enrichment  of  lipogenic  [^H]-NADPH  by  [3-^H]glucose  estimated  by  a  model  for  saturated  fatty  acid  synthesis  (ISA)  in  parental  H1299  and  A549  cells  following 
incubation  with  tracer  for  72  hr. 

Data  plotted  in  (A)-(D)  represent  mean  ±  SD  of  at  least  three  biological  replicates.  Eor  (E),  data  presented  are  mean  ±  95%  confidence  interval  of  at  least  three 
biological  replicates. 


the  reverse  cycle  may  be  used  to  produce  mitochondrial 
NADPH.  Metabolic  cycles  such  as  this  utilize  compartment-spe¬ 
cific  enzymes,  and  existing  methods  fortracing  metabolism  rely 
on  breaking  apart  cells  and  pooling  metabolites  from  all  com¬ 
partments,  making  it  impossible  to  reliably  distinguish  the  net  re¬ 
action  flux  through  each  enzyme  or  pathway. 

RESULTS 

Tracing  NADPH  with  ^H-Labeled  Glucose 

Because  reaction  mechanisms  involving  pyridine  nucleotides 
transfer  electrons  as  a  hydride  (H“)  ion,  isotope-labeled 
hydrogen  atoms  can  be  used  to  follow  electron  movement  in 
these  reactions  (Katz  et  al.,  1965;  Rendina  et  al.,  1984).  The 
transfer  of  and  can  also  be  used  to  observe  redox  reac¬ 
tions  in  central  carbon  metabolism,  an  approach  that  has  been 
used  to  generate  insight  into  NAD(P)H  metabolism  in  eukaryotic 
cells  (Ben-Yoseph  et  al.,  1994;  Ruhl  et  al.,  2012).  Glucose  is  the 
primary  carbon  source  for  glycolysis  and  the  oxidative  PPP  in 
mammalian  cells,  with  the  latter  pathway  representing  an  impor¬ 
tant  source  of  cytosolic  NADPH.  Nonlabile  hydrogen  atoms  on 
specific  glucose  carbons  (the  1  and  3  positions,  respectively) 
are  transferred  to  NADPH  by  the  oxidative  PPP  enzymes 


glucose-6-phosphate  dehydrogenase  (G6PD)  and  6-phospho- 
gluconate  dehydrogenase  (6PGD).  The  hydrogen  atom  on  car- 
bon-3  of  glucose  (which  becomes  carbon-1  of  dihydroxyace- 
tone  phosphate  in  glycolysis)  can  exchange  with  water  during 
isomerization  to  glyceraldehyde-3-phosphate  (GAP)  by  those 
phosphate  isomerase  (TPI)  (Katz  et  al.,  1 965).  This  prevents  con¬ 
founding  labeling  of  downstream  metabolites  including  TCA  cy¬ 
cle  intermediates,  suggesting  that  tracing  this  hydrogen  atom 
could  provide  a  means  of  quantifying  the  contribution  of  6PGD 
to  the  cellular  NADPH  pool  (Figure  1A)  (Katz  et  al.,  1966;  Katz 
and  Rognstad,  1978).  To  test  this  possibility,  we  cultured 
HI  299  non-small-cell  lung  cancer  cells  in  the  presence  of 
[3-^H]glucose  and  observed  labeling  of  NADPH  using  LC/MS- 
MS  (Figure  IB).  The  rapid  turnover  of  NADPH  allows  labeling 
from  [3-^H]glucose  to  reach  isotopic  steady  state  within 
30  min,  as  evidenced  by  the  lack  of  any  increased  label  incorpo¬ 
rated  into  NADPH  after  culturing  cells  in  the  presence  of  [3-^H] 
glucose  for  24  hr  (Figure  1 B).  NADPH  has  two  hydrogens  that 
can  be  transferred  when  it  acts  as  an  electron  donor.  Once 
labeled,  either  the  labeled  or  unlabeled  hydrogen  atom  can  be 
transferred  depending  on  the  stereospecificity  of  downstream 
NADPH-utilizing  enzymes  (You,  1985).  Transfer  of  the  unlabeled 
hydride  from  labeled  NADPH  generates  labeled  NADP"^  (see 
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Figure  S1 A  available  online),  and  subsequent  labeling  of  the  sec¬ 
ond  hydrogen  on  NADP"^  yields  NADPH  heavy  by  two  mass  units 
(M2)  at  later  time  points  (Figure  1 B).  Some  labeling  of  ribose-5- 
phosphate  and  ribulose-5-phosphate  was  also  observed  at 
late  time  points,  presumably  via  flux  through  the  nonoxidative 
PPP  (Figure  S1 B).  This  label  incorporation  into  the  ribose  moi¬ 
eties  of  NADP(H)  could  account  for  a  minor  portion  of  the  isotope 
enrichment  observed  at  24  hr,  as  suggested  by  the  small  amount 
of  M2  labeling  of  NADP"^  (Figure  SI  A),  but  these  atoms  would 
not  be  subject  to  transfer  in  downstream  reactions  utilizing 
NADPH  as  a  cofactor.  [3-^H]glucose  does  not  contribute  to 
NAD(H)  except  by  incorporation  into  the  ribose  moieties,  and 
oniy  a  small  amount  of  NAD(H)  is  labeled  at  late  time  points  (Fig¬ 
ure  SI  C),  arguing  that  the  majority  of  NADPH  labeling  from  [3-^H] 
glucose  reflects  hydride  transfer.  Importantly,  metabolites  in 
lower  glycolysis  such  as  pyruvate  and  lactate  are  not  labeled 
from  [3-^H]glucose  (see  percent  values  in  Figure  1A),  implying 
that  the  presence  of  label  on  downstream  metabolites  must  arise 
as  a  result  of  hydride  transfer  from  labeled  NADPH.  [1-^H] 
glucose  labels  NADPH  via  G6PD  in  the  oxidative  PPP  (Figures 
SI  D  and  SI  E),  and  similar  to  [3-^H]glucose,  this  label  is  detected 
on  NADP"^  and  PPP  intermediates  (Figures  S1E  and  S1F).  How¬ 
ever,  deuterium  present  on  carbon-1  of  [1  -^Hjglucose  can  be  lost 
due  to  reversibility  of  phosphoglucose  isomerase,  resulting  in 
less  glucose-6-phosphate  labeling  from  [1-^H]glucose  in  cells 
compared  to  labeling  from  [3-^H]glucose  (Figure  S1G)  (Ben-Yo- 
seph  et  al.,  1994;  Hellerstein  et  al.,  1986;  Katz  and  Rognstad, 
1976).  In  contrast  to  [3-^H]glucose,  the  deuterium  isotope  from 
carbon-1  of  glucose  does  not  exchange  with  water  in  any  of 
the  reactions  of  glycolysis  and  is  retained  on  carbon  entering 
the  TCA  cycle  (Figure  SID).  Therefore,  deuterium  from  [1-^H] 
glucose  has  the  potential  to  label  downstream  metabolites  by 
either  hydride  transfer  from  NADPH  or  label  retention  on  the  car¬ 
bon  from  glucose.  As  a  result,  use  of  [3-^H]glucose  is  preferable 
for  tracing  NADPH  produced  by  the  oxidative  PPP. 

The  NADPH  labeling  we  observe  is  a  subset  of  the  total  cellular 
pool.  To  gain  insights  into  compartment-specific  redox  reac¬ 
tions,  we  next  quantified  ^H  enrichment  in  specific  metabolites 
downstream  of  NADPH-dependent  reactions.  For  example,  fatty 
acid  and  cholesterol  synthesis  occur  specifically  in  the  cytosol 
and  require  NADPH.  When  HI 299  cells  are  cultured  with  [3-^H] 
glucose  for  72  hr  to  allow  for  accumulation  of  new  lipid  mole¬ 
cules,  we  detected  significant  labei  on  newly  synthesized  fatty 
acids,  including  myristate  (Cl  4:0),  paimitate  (Cl  6:0),  and  stea¬ 
rate  (Cl 8:0),  as  well  as  cholesterol  (Figures  1C  and  ID).  Impor¬ 
tantly,  no  label  from  [3-^H]glucose  was  detected  on  citrate  (see 
percent  value  in  Figure  1A),  suggesting  that  isotope  enrichment 
on  lipids  was  from  the  NADPH  pool  that  was  labeled  by  the 
oxidative  PPP.  We  also  detected  labeling  of  fatty  acids  from 
[1-^H]-glucose  (Figure  S1H);  however,  tracing  of  G6PD-derived 
NADPH  is  complicated  by  deuterium  from  [1-^H]glucose  being 
retained  on  citrate  and  lipogenic  acetyl-CoA  (Figure  SI  D).  We 
next  applied  isotopomer  spectral  analysis  (ISA)  to  estimate  the 
contribution  of  [3-^H]glucose  to  the  lipogenic  NADPH  pool 
used  for  paimitate  synthesis,  as  14  NADPH  molecules  are 
required  during  the  production  of  one  paimitate  molecule  (Khar- 
roubi  et  al.,  1992;  Metallo  et  al.,  2012).  The  ISA  model  includes 
two  parameters  representing  the  deuterium  enrichment  of  the 


NADPH  pool  and  the  percentage  of  paimitate  that  was  synthe¬ 
sized  de  novo  (Figure  S2).  Using  this  method  we  estimated  the 
enrichment  of  lipogenic  NADPH  from  [3-^H]glucose  ranged 
from  12%-20%  in  A549  and  H1299  cells  (Figure  IE).  ISA 
modeling  of  other  saturated  fatty  acids  (e.g.,  myristate  and  stea¬ 
rate)  yielded  similar  estimations  for  the  enrichment  of  lipogenic 
NADPH  from  [3-^H]glucose  (Figure  IE). 

Use  of  Glucose  to  Trace  NADH  Metabolism 

To  maintain  flux  through  glycolysis,  cytosolic  NAD"^  pools  are  re¬ 
generated  primarily  by  three  enzymes:  lactate  dehydrogenase 
(LDH),  malate  dehydrogenase  (MDH),  and/or  the  glycerol  phos¬ 
phate  shuttle  (Glyc3PDH)  (Lunt  and  Vander  Heiden,  201 1 ;  Met¬ 
allo  and  Vander  Heiden,  2013).  Distinct  hydrogen  atoms  on 
glucose  are  transferred  to  NAD"^  during  glycolysis  via  glyceralde- 
hyde  phosphate  dehydrogenase  (GAPDH).  In  theory,  up  to  half  of 
the  hydrogen  transferred  to  NADH  via  GAPDH  comes  from  car¬ 
bon  four  of  glucose;  however,  exchange  with  water  in  the 
aldolase  and  TPI  reactions  decreases  the  net  contribution  of 
this  hydrogen  atom  to  NADH  (Go  et  al.,  2009)  (Figure  2A). 
Upon  culturing  A549  and  HI 299  cells  with  [4-^H]glucose,  signif¬ 
icant  labeling  of  lactate,  malate,  and  glycerol  3-phosphate  was 
observed  (Figure  2B).  Label  was  detected  on  GAP  in  A549  cells; 
however,  the  level  of  GAP  in  HI  299  cells  was  below  the  limit  of 
detection.  In  addition,  no  label  was  detected  on  metabolites  in 
lower  glycolysis  including  PEP,  3PG,  and  pyruvate.  This  pattern 
fits  with  known  reactions  using  NADH  in  central  carbon  meta¬ 
bolism  (Figure  2A)  and  suggests  that  [4-^H]glucose  can  be 
used  to  label  the  NADH  pool  produced  by  glycolysis  in  cells. 
Consistent  with  these  findings,  we  observed  rapid  labeling 
from  [4-^H]glucose  on  NADH  (Figure  2C),  as  well  as  label  on 
NAD+  (Figure  S3A)  arising  from  donation  of  the  unlabeled  hy¬ 
dride  from  Ml  labeled  NADH. 

Interestingly,  we  observed  isotope  incorporation  into  fatty  acid 
pools  from  [4-^H]glucose  (Figure  2D),  suggesting  that  some  label 
from  NADH  transfers  to  cytosolic  NADPH  through  an  unknown 
mechanism  (Figure  S3F).  While  deuterium  label  from  [4-^H] 
glucose  was  detected  on  aspartate,  citrate,  and  isocitrate  due 
to  the  symmetry  of  fumarate  (Figures  S3B-S3D),  carbons 
labeled  in  this  manner  do  not  contribute  to  lipogenic  acetyl- 
CoA,  demonstrating  that  the  observed  fatty  acid  labeling  is 
derived  from  hydride  transferring  from  NAD(H)  to  cytosolic 
NADP(H).  We  also  observed  some  isotope  enrichment  on 
NADP"^  and  NADPH  from  [4-^H]glucose  (Figures  2E  and  S3A); 
however,  ribose  5-phosphate  and  ribulose  5-phosphate  are 
also  labeled  from  [4-^H]glucose  (Figure  S3E).  Furthermore,  these 
direct  measurements  of  total  cellular  NADP(H)  cannot  distin¬ 
guish  between  cytosolic  and  mitochondrial  pools,  highlighting 
the  need  for  methods  to  elucidate  compartment-specific 
NADP(H)  pools. 

A  Reporter  System  to  Trace  Compartmentalized 
Sources  of  NADPH 

The  above  data  demonstrate  that  we  can  observe  cytosolic 
production  of  NADPH  and  NADH  in  intact  cells.  Although 
we  were  able  to  quantify  the  contribution  of  oxidative  PPP 
enzymes  to  the  lipogenic  NADPH  pooi,  deuterium  tracing 
alone  cannot  distinguish  other  compartmentalized  sources  of 
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Figure  2.  Use  of  Glucose  to  Label  NADH 

(A)  Atom-transition  map  depicting  a  model  of  deuterium  transfer  from  [4-^H]glucose  through  glycolysis  and  NAD+-dependent  shuttle  systems  (malate  dehy¬ 
drogenase  MDH,  glycerol  3-phosphate  dehydrogenase  GlySPDH,  and  lactate  dehydrogenase  LDH).  Open  large  circles  represent  carbon,  and  small  green  circles 
indicate  deuterium  label  from  [4-^H]glucose. 

(B)  Labeling  of  glycolytic  intermediates  from  [4-^H]glucose  in  A549  (left  panel)  and  HI  299  (right  panel)  cells.  Glyceraldehyde-3-phosphate  (GAP)  was  below  the 
limit  of  detection  in  HI  299  cells,  indicated  by  asterisk. 

(C)  Labeling  of  NADH  from  [4-^H]glucose  in  parental  HI 299  cells  overtime. 

(D)  Palmitate  labeling  from  [4-^H]glucose  in  A549  and  HI  299  cells  following  incubation  with  tracer  for  72  hr. 

(E)  Time  course  labeling  of  NADPH  from  [4-^H]glucose  in  parental  H1299.  Data  presented  are  mean  ±  SD  of  at  least  three  biological  replicates. 


NADPH.  Therefore,  we  sought  to  develop  a  reporter  system 
that  can  detect  pathway-specific  NADPH  production  in  dif¬ 
ferent  subcellular  compartments.  To  accomplish  this,  we  took 
advantage  of  the  neomorphic  mutant  IDH  enzymes  that  pro¬ 
duce  (D)2-hydroxyglutarate  (2HG)  from  aKG.  This  reaction 
reduces  aKG  by  transferring  a  hydride  from  NADPH  to  form 
2HG.  As  2HG  is  a  xenometabolite  that  is  only  present  at  very 
low  levels  in  most  cells  (Matsunaga  et  al.,  2012),  it  can  be 
used  as  an  endproduct  readout.  By  applying  specific  metabolic 
^H-tracers  to  cells  and  measuring  enrichment  of  2HG  produced 
by  ectopically  expressed  mutant  IDH1  (cytosol)  or  IDH2  (mito¬ 
chondria),  we  reasoned  that  pathway-specific  information  on 
NADPH  metabolism  in  each  compartment  could  be  obtained 
(Figure  3A). 

We  generated  HI  299  and  A549  cell  lines  that  express  epitope- 
tagged  mutant  IDH1-R132H  (mtIDHI-C)  or  mutant  IDH2-R172K 
(mtlDH2-M)  in  a  doxycycline-dependent  manner  (Figure  S4A). 
We  employed  a  weak  promoter  to  minimize  effects  on  endo¬ 
genous  IDH  metabolism.  Indeed,  FLAG-tagged  mtIDHI-C  was 
expressed  at  levels  that  were  not  detectable  with  an  antibody 


recognizing  wild-type  IDH1  enzyme  in  these  cells  (Figure  S4A). 
mtIDHI-C  is  expected  to  be  expressed  in  the  cytoplasm  and 
mtlDH2-IVI  in  the  mitochondria,  and  we  confirmed  the  localiza¬ 
tion  of  each  using  cell  fractionation  and  western  blotting  (Fig¬ 
ure  3B).  We  also  confirmed  that  the  FLAG-tagged  mutant  IDH 
enzymes  produce  2HG  in  a  doxycycline-dependent  manner  in 
both  HI 299  and  A549  cell  lines  (Figure  3C).  Interestingly, 
mtIDHI-C  produced  less  2HG  than  mtlDH2-M  in  cells,  consis¬ 
tent  with  observations  that  ectopically  expressed  IDH2  mutants 
produce  more  2HG  than  IDH1  mutants  due  to  their  mitochondrial 
localization  (Ward  et  al.,  2013).  In  all  cases,  2HG  levels  were  far 
below  those  observed  in  tumor  cell  lines  expressing  endogenous 
mutations  in  IDH1  (R132C/-I-,  HT1080)  or  IDH2  (R172S/-r, 
SW1353)  (Figure  3D).  Introduction  of  mutant  IDH  enzymes  could 
impact  NADPH  or  TCA  metabolism;  however,  2HG  production 
flux  observed  in  cell  lines  expressing  IDH1  mutants  at  higher 
levels  is  small  relative  to  other  aKG-dependent  reaction  fluxes, 
suggesting  mutant  IDH  expression  has  a  minimal  direct  impact 
on  aKG  pools  (Grassian  et  al.,  2014).  Furthermore,  no  significant 
change  in  [3-^H]glucose  contribution  to  cytosolic  NADPH  was 
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Figure  3.  Generation  and  Characterization  of  Cell  Lines  Expressing  Inducible  Mutant  IDH 

(A)  Schematic  demonstrating  the  transfer  of  deuterium  (red  dots)  from  NADPH  to  2HG  via  the  reaction  cataiyzed  by  mutant  iDH  enzymes.  mtiDHI  is  locaiized  to 
the  cytopiasm,  and  mtiDH2  is  locaiized  to  the  mitochondria.  By  expressing  compartment-specific  mutant  enzymes  and  combining  this  with  deuterated  giucose 
tracer,  it  is  possible  to  track  2HG  production  and  therefore  the  source  of  the  NADPH  used  to  make  the  2HG  in  the  cytosoi  or  the  mitochondria. 

(B)  Mutant  iDH1-R132H  is  localized  to  the  cytosol  (mtIDHI-C),  and  mutant  IDH2-R172K  is  localized  to  the  mitochondria  (mtlDH2-M)  in  H1299  and  A549  cells. 
Cells  were  transduced  with  lentiviral  constructs  containing  cDNA  encoding  C-terminal  FLAG-tagged  IDH1-R132H  or  IDH2-R172K  under  the  control  of  a 
doxycycline  (dox)-inducible  promoter.  Once  stable  cell  lines  were  established,  cells  were  treated  with  0.1  ng/mL  dox  for  24  hr.  Protein  expression  was  analyzed 
by  cellular  fractionation  and  western  blotting  using  antibodies  against  FLAG,  cytochrome  c  (mitochondrial-specific  marker),  and  Hsp70  (cytoplasmic-specific 
marker).  C,  Supernatant-100  fraction  (cytoplasm);  M,  mitochondria.  Samples  from  mtiDHI -C  and  mtlDH2-M  cells  were  analyzed  on  separate  gels. 

(C)  Cell  lines  expressing  inducible  IDH  mutants  produce  2HG  in  a  dox-dependent  manner.  H1299  and  A549  cells  stably  expressing  inducible  mtIDHI-C  or 
mtlDH2-M  constructs  were  treated  with  doxycycline  (0.1  ng/mL)  for  24  hr.  Amounts  of  2HG  (total  ion  counts,  TIC)  are  shown  relative  to  GFP  control  cells  treated 
with  vehicle. 

(D)  2HG  production,  as  measured  by  2HG/aKG  ratio,  is  much  higher  in  cell  lines  harboring  endogenous  mutations  in  IDH1  (R1 32C/-t-,  HT 1 080)  and  IDH2  (R1 72S/-I-, 
SW1353)  than  cells  expressing  mtiDHI -C  and  mtlDH2-M. 

(E)  NADPH  produced  by  the  pentose  phosphate  pathway  (6PGD)  is  cytosolic.  Cells  were  cultured  in  [3-^H]-glucose  (10  mM)  for  24  hr  before  adding  dox 
(0.1  rig/mL)  for  24  hr  to  induce  mutant  IDH  expression,  and  amount  of  Ml  label  (%)  from  [3-^H]-glucose  incorporated  into  2HG  and  aKG  was  measured. 

(F)  NADH  supports  NADPH  production  in  the  mitochondria.  Cells  were  incubated  with  10  mM  [4-^H]-glucose  for  24  hr  and  treated  and  analyzed  as  in  E.  Data 
represent  mean  ±  SEM  of  at  least  three  biological  replicates. 


observed  in  A549  cells  following  mtiDHI  -C  or  mtlDH2-M  expres¬ 
sion  (Figure  S4G). 

Although  doxycycline  can  affect  the  metabolism  and  prolifer¬ 
ation  of  some  mammalian  cancer  cell  lines  in  culture  (Abler 
et  al.,  2013),  we  saw  no  doxycycline-dependent  changes  in  the 


abundance  of  central  carbon  metabolites  (Figure  S4B)  or  in  the 
proliferation  rate  (Figure  S4C)  of  A549  or  HI 299  cells.  Impor¬ 
tantly,  cells  expressing  doxycycline-inducible  GFP  also  showed 
no  changes  in  metabolite  pool  sizes  or  proliferation  rates  (Figures 
S4B  and  S4C),  indicating  that  when  added  at  this  concentration 
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doxycycline  does  not  significantly  affect  metabolism  in  this  sys¬ 
tem.  In  addition,  we  observed  no  significant  differences  in  pool 
sizes  of  NAD-",  NADH,  NADP-",  or  NADPH  in  H1299  mtlDH1-C 
and  mtlDH2-M  cells  following  the  addition  of  doxycycline  for 
24  hr,  suggesting  that  doxycycline-dependent  production  of 
2HG  was  not  altering  the  availability  of  these  cofactors  for  use 
in  other  redox  reactions  (Figure  S4D).  This  is  supported  by  the 
reported  kcat  of  mutant  IDH1  enzymes  being  small  relative  to 
wild-type  IDH1  (Dang  et  al.,  2009)  and  suggests  that  any  direct 
effects  of  these  enzymes  on  cellular  redox  state  are  minimal. 

Validation  of  Compartment-Specific  Cofactor  Tracing 

To  validate  the  ability  of  this  system  to  trace  compartment-spe¬ 
cific  NADPH  metabolism,  we  induced  expression  of  the  mutant 
IDH  enzymes  in  cells  cultured  in  the  presence  of  [3-^H]glucose 
and  measured  enrichment  of  in  the  2HG  pool.  In  order  to 
ensure  that  the  cells  were  at  or  near  isotopic  steady  state  prior 
to  induction  of  mutant  IDH  expression,  the  cells  were  incubated 
with  tracer  for  24  hr  prior  to  the  addition  of  doxycycline.  Consis¬ 
tent  with  [3-^H]glucose  producing  cytosolic  NADPH  via  the 
oxidative  PPP,  2HG  was  only  significantly  labeled  from  [3-^H] 
glucose  in  the  mtlDH1-C  cell  lines  and  not  in  the  mtlDH2-M 
cell  lines  (Figure  3E).  Importantly,  little  to  no  label  was  observed 
on  aKG  under  these  conditions,  ensuring  that  label  on  2HG  was 
a  direct  result  of  hydride  ion  transfer  from  NADPH  by  the  mutant 
enzyme.  We  next  asked  whether  label  from  [4-^H]glucose  was 
incorporated  into  2HG  by  either  mtlDH1-C  or  mtlDH2-M. 
Notably,  more  2HG  was  labeled  from  [4-^H]glucose  in  mtlDH2- 
M  cells  than  in  mtlDH1-C  cells  (Figure  3F).  These  data  suggest 
that  transfer  of  H“  from  NADH  to  NADPH  occurs  through  a 
mitochondrial  intermediate  (e.g.,  malate)  or  via  nicotinamide 
nucleotide  transhydrogenase  (NNT),  and  transfer  of  reducing 
equivalents  from  NADH  to  NADPH  mostly  supports  the  mito¬ 
chondrial  NADPH  pool.  Similar  results  were  observed  in  cell  lines 
with  endogenous,  heterozygous  IDH1  and  IDH2  mutations  (Fig¬ 
ures  S4E  and  S4F). 

In  vitro  steady-state  enzyme  kinetics  experiments  have 
demonstrated  that  rate  constants  for  reactions  involving  ^H 
transfer  can  be  lower  compared  to  studies  conduoted  with  unla¬ 
beled  substrate  (Rendina  et  al.,  1984).  This  phenomenon  has 
been  an  invaluable  tool  for  elucidating  biochemical  reaction 
mechanisms,  including  many  of  the  reactions  responsible  for 
label  transfer  in  our  system.  However,  given  the  diverse  means 
through  which  metabolism  is  regulated  and  the  ohallenges 
associated  with  understanding  which  enzymatic  steps  are  rate 
limiting  for  pathways  in  intact  cells,  the  relevance  of  isotope  ef¬ 
fects  to  intracellular  metabolic  fluxes  is  not  clear.  To  determine 
the  significance  of  this  phenomenon  in  our  system,  we  cultured 
cells  in  different  ratios  of  [3-^H]gluoose  and  unlabeled  glucose 
and  measured  downstream  labeling  of  lipogenio  NADPH  and 
2HG  in  HI 299  mtIDHI-C  oells  (Figure  4A).  Similarly,  we  titrated 
[4-^H]glucose  with  unlabeled  substrate  in  HI  299  mtlDH2-M 
cells  and  observed  whether  label  transfer  to  lactate,  malate, 
aspartate,  fumarate,  citrate,  and  2HG  was  affected  by  different 
amounts  of  labeled  substrate  (Figure  4B).  We  reasoned  that  if 
reaction  rates  are  affected  by  the  presence  of  ^H  isotopes,  the 
use  of  unlabeled  substrates  would  be  favored  and  less  relative 
transfer  of  label  would  be  observed  as  unlabeled  substrate  is 


titrated  into  the  medium.  However,  in  all  cases,  transfer  of  label 
from  either  tracer  decreased  linearly  as  the  tracer  was  diluted 
with  unlabeled  substrate,  suggesting  that  kinetic  isotope  effeots 
minimally  impact  the  results  of  these  experiments  (Figures  4A 
and  4B).  A  linear  decrease  in  lipogenic  NADPH  and  2HG  labeling 
was  also  observed  in  endogenous  mtIDHI  cells  (HT1080)  that 
exhibit  much  higher  rates  of  2HG  production  when  we  titrated 
[3-^H]glucose,  further  supporting  the  notion  that  isotope  effects 
minimally  affect  substrate  fluxes  through  the  reactions  we  traced 
in  intact  cells  (Figure  4C). 

Characterizing  Serine/Glycine  Metabolism  in  the 
Cytosol  and  Mitochondria 

We  next  sought  to  use  our  reporter  system  to  examine  a 
compartmentalized  metabolic  cycle.  Reactions  that  make  up 
folate-mediated  one-carbon  metabolism  exist  in  both  the 
cytosol  and  the  mitochondria,  although  it  is  not  clear  from  current 
literature  whether  these  are  linked  in  a  cycle  and/or  in  which 
direction  the  reactions  prooeed  (Anderson  et  al.,  2011;  Nilsson 
et  al.,  2014;  Tedeschi  et  al.,  2013;  Tibbetts  and  Appling,  2010). 
Serine  and  glyoine  interconversion  via  serine  hydroxymethyl- 
transferase  (SHMT)  has  been  observed  in  cultured  cells  (Jain 
et  al.,  2012;  Levintow  and  Eagle,  1961;  Perry  et  al.,  2007),  but 
tracing  is  unable  to  ascertain  the  directionality,  compartmen- 
talization,  and  interconneotivity  of  this  prooess.  Indeed,  upon 
culture  with  [U-^^Calserine,  we  observed  significant  interconver¬ 
sion  of  serine  and  glyoine  in  A549  mtIDHI-C  and  mtlDH2-M  cells 
(Figure  5A).  The  reactions  catalyzed  by  MTHFD1  (cytosolic)  and 
MTHFD2/MTHFD2L  (mitochondrial)  utilize  NAD(P)H  (Figure  5B); 
therefore,  we  hypothesized  that  ^H  serine  and  glycine  tracing 
in  combination  with  our  oompartment  reporter  system  would 
enable  us  to  experimentally  determine  the  direotion  of  serine- 
glycine  exchange  reactions  in  the  cytosol  and  mitoohondria. 
Hydrogens  on  carbon-3  of  serine  are  transferred  to  5,10-methyl- 
enetetrahydrofolate  (5,10-methyleneTHF)  and  subsequently  to 
NADPH  via  MTHFD1/2  (Figure  5B).  Additionally,  the  glycine 
cleavage  system  (GCS)  exists  in  the  mitochondria  and  could 
transfer  hydrogen  from  carbon  two  of  glycine  to  5,10-methy- 
lene-THF  and  generate  NADPH  (Figure  5B)  (Kikuchi  et  al.,  2008). 

To  study  these  compartment-specific  pathways,  we  cultured 
A549  mtIDHI-C  and  mtlDH2-M  cells  with  either  [3,3-^H2]serine 
or  [2,3,3-^H3]serine  and  unlabeled  glycine  or  [2,2-^H2]glycine 
and  unlabeled  serine  and  measured  incorporation  of  ^H  in  cyto- 
solio  or  mitochondrial  2HG,  respeotively.  Strikingly,  we  detected 
label  from  [3,3-^H2]serine  and  [2,3,3-^H3]serine  on  2HG  only 
in  mtlDH2-M  cells,  strongly  suggesting  that  serine  to  glycine 
conversion  occurs  primarily  in  the  mitochondria  in  these  cells 
with  the  MTHFD2/MTHFD2L  reaction  operating  oxidatively  (Fig¬ 
ures  5C  and  S5A-S5C).  We  did  not  observe  labeling  of  2HG  from 
[2,2-^H2]glycine  in  oells  expressing  either  mutant  IDH  (Figure  5C), 
indicating  that  either  the  majority  of  mitochondrial  glyoine  is 
generated  by  SHIVIT2  (rather  than  glycine  import)  or  the  label  is 
lost  in  the  GCS.  Consistent  with  the  lack  of  label  transfer  from 
^H-labeled  serine  or  glycine  to  2HG  in  mtIDHI-C  cells,  we  de¬ 
tected  minimal  contribution  of  these  tracers  in  the  lipogenic 
NADPH  pool  (Figure  5D).  To  further  oonfirm  the  direction  of 
MTHFD1,  we  cultured  A549  mtIDHI-C  and  mtlDH2-M  cells 
with  [3-^H]glucose,  which  specifically  labels  cytosolic  NADPH. 
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[4-^H]glucose  media  enrichment  (%) 


Figure  4.  Kinetic  Isotope  Effect  Minimally  Affects  [3-^H]Glucose  and  [4-^H]Glucose  Metabolism 

(A)  [3-^H]glucose  was  titrated  with  uniabeied  glucose  and  added  to  HI  299  ceiis  expressing  mtiDHI  -C.  Contribution  from  [3-^H]giucose  to  iipogenic  NADPH  (left 
panei,  normalized  to  contribution  at  100%  [3-^H]giucose  media  enrichment)  and  2HG  (right  panel)  was  measured.  Dashed  lined  (left  panel)  represents  1:1 
contribution  of  [3-^H]glucose  to  Iipogenic  NADPH  to  enrichment  of  [3-^H]glucose  in  media. 

(B)  [4-^H]glucose  was  titrated  with  unlabeled  glucose  and  added  to  H1299  cells  expressing  mtlDH2-M.  Labeling  from  [4-^H]glucose  on  lactate  and  malate 
(left  panel);  aspartate,  citrate,  and  fumarate  (middle  panel);  and  2HG  (right  panel)  was  measured. 

(C)  [3-^H]glucose  was  titrated  with  unlabeled  glucose  in  HT1080  cells  harboring  endogenous  IDH1  mutations  (R132C/+).  Contribution  from  [3-^H]glucose 
to  Iipogenic  NADPH  (left  panel,  normalized  to  contribution  at  100%  [3-^H]glucose  media  enrichment)  and  Ml  labeling  on  2HG  (right  panel)  was  quantified  at 
0%,  25%,  50%,  75%,  and  100%  dilution  with  unlabeled  glucose  (left  panel)  in  HT1080  cells  cultured  with  [3-^H]glucose  diluted  with  unlabeled  glucose  at 
0%,  25%,  50%,  75%,  and  100%  enrichment. 

Data  presented  are  shown  as  mean  ±  SD  of  three  biological  replicates  for  (A)  (right  panel),  (B),  and  (C)  (right  panel).  Data  represent  mean  ±  95%  confidence  interval 
for  at  least  three  biological  replicates  for  (A)  (left  panel)  and  (C)  (left  panel). 


We  observed  transfer  of  H  from  [3-  H]glucose  onto  serine,  sug¬ 
gesting  that  cyfosoiic  MTHFD1  can  operate  in  the  reductive 
direction  in  these  ceiis  (Figure  5E).  The  iack  of  giycine  iabeling 
from  [3-^H]giucose  confirms  iabei  transfer  to  serine  was  ob¬ 
tained  from  5,10-methyieneTHF  (Figures  5E  and  S5D).  The  iack 
of  iabeiing  from  either  [2,3,3-^Fi3]serine  or  [3,3-^H2]serine  on 
either  fatty  acids  or  2FiG  produced  by  mtIDFII  -C  is  aiso  consis¬ 
tent  with  minimai  contribution  of  MTHFD1  to  the  cytopiasmic 
NADPH  pooi  in  these  ceiis,  aithough  channeiing  to  other  reac¬ 
tions  cannot  be  ruied  out  by  these  methods.  Shuttiing  of  serine 
iabeied  by  [3-^H]giucose  by  MTHFD1/SHMT1  into  the  mito¬ 
chondria  for  cataboiism  by  MTHFD2/SHMT2  may  account  for 
the  smaii  amount  of  iabei  (<1  %)  we  observe  on  mitochondriai 
2HG  in  A549  mtiDH2-IVl  ceiis  cuitured  with  [3-^H]giucose  (Fig¬ 
ure  3E).  Coiiectiveiy,  these  data  provide  direct  evidence  that 
serine  metaboiism  can  contribute  to  regenerating  mitochondriai 
NADPH  in  ceiis. 

DISCUSSION 

We  have  deveioped  a  system  that  can  distinguish  compartmen¬ 
taiized  poois  of  NADPH,  demonstrating  the  directionality  and 
interconnectivity  of  serine/glycine  metabolism  in  the  cytosol 
and  mitochondria  of  intact  cells.  In  the  cells  studied,  conversion 
of  serine  to  glycine  occurs  primarily  in  the  mitochondria  with 


the  reaction  catalyzed  by  MTHFD2/IVITHFD2L  contributing  to 
NADPH  production  in  this  compartment.  Interestingly,  label 
transfer  from  both  [2,3,3-^H3]serine  and  [3,3-^H2]serine  to  mito¬ 
chondrial  2HG  was  observed,  suggesting  that  serine  meta¬ 
bolism  by  SHMT2  is  a  contributor  to  the  mitochondrial  NADPH 
and  glycine  pools.  These  data  also  provide  direct  experimental 
support  for  the  hypothesis  that  the  cytoplasmic  source  of 
formate  used  for  purine  synthesis  can  be  mitochondrially 
derived  in  some  cells.  Previous  efforts  to  ascertain  directionality 
of  folate-mediated  one-carbon  metabolism  have  been  unable  to 
distinguish  between  compartments,  relying  on  expression  data 
(Nilsson  et  al.,  2014),  mathematical  modeling  (Scotti  et  al., 
2013;  Tedeschi  et  al.,  2013),  or  isolated  mitochondrial  prepara¬ 
tions  (Barlowe  and  Appling,  1988).  The  importance  of  distin¬ 
guishing  compartmentalized  redox  pathways  is  highlighted  by 
the  large  number  of  potential  pathways  that  have  been  impli¬ 
cated  in  the  shuttling  of  reducing  equivalents  between  the 
cytosol  and  mitochondria.  For  instance,  compartment-specific 
metabolic  cycling  through  citrate/aKG  (Sazanov  and  Jackson, 
1994;  Ward  et  al.,  2010),  malate/pyruvate  (Jiang  et  al.,  2013; 
Son  et  al.,  2013),  proline  (Hagedorn  and  Phang,  1983;  Nilsson 
et  al.,  2014),  and  serine  (Tibbetts  and  Appling,  2010)  have 
been  suggested  to  be  important  for  mammalian  cell  physiology. 
Although  carbon  tracing  is  increasingly  combined  with  genetic 
approaches  to  implicate  a  role  for  compartment-specific 
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Figure  5.  Characterizing  Serine/Glycine  Metabolism  in  the  Cytosol/Mitochondria 

(A)  Serine  (left  panel)  and  glycine  (right  panel)  labeling  in  A549  mtIDHI-C  and  mtlDH2-M  cells  cultured  with  [U-^^Calserine.  Cells  were  incubated  with  [U-^^Ca] 
serine  for  24  hr  prior  to  dox-induction  (0.1  rig/mL)  for  an  additional  48  hr.  Middle  panel  demonstrates  interconversion  of  serine  and  glycine  by  SHMT.  Data  plotted 
represent  mean  ±  SD  for  three  biological  replicates. 

(B)  A  schematic  of  folate-mediated  one-carbon  metabolism  in  cytosolic  and  mitochondrial  compartments  catalyzed  via  SHMT  and  MTHFD.  Deuterium  transfer 
from  [3,3-^H2]serine  is  shown  for  pathways  containing  SHMT  and  MTHFD  and  is  indicated  by  small  red  or  blue  circles  for  cytosolic  and  mitochondrial  isozymes, 
respectively.  The  extra  deuterium  on  [2,3,3-^H3]serine  is  indicated  by  an  orange  (cytosolic)  or  a  turquoise  (mitochondrial)  circle.  Deuterium  transfer  from  [2,2-^H2] 
glycine  is  shown  for  the  glycine  cleavage  system  (GCS)  pathway  indicated  by  small  green  circles. 

(C)  2HG  labeling  from  [3,3-^H2]serine,  [2,3,3-^H3]serineor  [2,2-^H2]glycine  in  A549  mtlDH1-C  and  mtlDH2-M  cells.  Cells  were  incubated  with  either  tracer  for  24  hr 
prior  to  dox  induction  (0.1  ng/mL)for  an  additional  48  hr.  No  label  was  detected  on  2HG  in  mtIDHI-C  cells  from  either  [3,3-^H2]serineor  [2,3,3-^H3]serine,  nor  was 
label  detected  on  2HG  from  [2,2-^H2]glycine  in  mtIDHI-C  and  mtlDH2-M  cells  (indicated  by  asterisk). 

(D)  Fatty  acid  labeling  from  A549  mtlDHI-C  and  mtlDH2-M  cells  cultured  with  either  [3,3-^H2]serine,  [2,3,3-^H3]serine,  or  [2,2-^H2]glycine.  Cells  were  incubated 
with  tracer  for  24  hr  prior  to  dox  induction  (0.1  rrg/mL)  for  an  additional  48  hr. 

(E)  Serine  and  glycine  labeling  in  A549  mtIDHI-C  and  mtlDH2-M  cells  cultured  with  [3-^H]glucose.  Cells  were  incubated  with  tracer  for  24  hr  prior  to  dox  induction 
(0.1  pg/mL)  for  48  hr.  Data  represent  mean  ±  SEM  of  at  least  three  biological  replicates. 


isozymes  in  such  processes,  adaptation  to  genetic  depletion 
strategies  that  break  these  cycles  can  confound  interpretation. 
The  reporter  system  described  here  circumvents  these  issues 


for  reactions  involving  NADPH  by  providing  direct  visualization 
of  compartmentalized  reaction  activity  and  direction  in  intact 
cells. 
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other  subcellular  compartments  also  have  distinct  metabolic 
needs  in  eukaryotic  cells  and  could  be  probed  with  an  analogous 
approach  by  engineering  the  localization  of  mutant  IDH  enzymes 
and/or  monitoring  transfer  between  other  metabolites.  The 
endoplasmic  reticulum  (ER)  is  an  important  site  for  protein 
folding,  disulfide  bond  formation,  long  chain  fatty  acid  extension, 
and  sterol  reduction;  as  such,  the  NADP'^/NADPH  ratio  within  the 
ER  can  influence  diverse  cellular  and  physiological  processes 
(Banhegyi  et  al.,  2009;  Kardon  et  al.,  2008;  Szaraz  et  al.,  2010). 
This  approach  may  also  be  adapted  to  quantify  the  NADPV 
NADPH  ratio  in  particular  organelles  and  the  turnover  rate  of 
NADPH  in  specific  compartments.  Altering  the  cofactor  selec¬ 
tivity  of  the  mutant  IDH  enzyme,  or  relying  on  NADH-dependent 
production  of  another  xenometabolite  could  similarly  be  used  to 
visualize  compartmentalized  reactions  that  utilize  NADH.  Finally, 
these  data  may  be  integrated  with  compartmentalized  meta¬ 
bolic  flux  analysis  (MFA)  models  to  better  understand  the  role  of 
cofactor  metabolism  in  metabolic  engineering  applications  or 
disease  models  (Ruhl  et  al.,  2012).  Thus,  this  approach  opens 
up  new  avenues  to  observe  metabolic  processes  in  complex 
cells  and  improve  our  understanding  of  metabolism  in  normal 
and  disease  states. 

EXPERIMENTAL  PROCEDURES 
Cell  Culture  and  Isotopic  Labeling 

All  cell  lines  were  maintained  in  DMEM  supplemented  with  10%  FBS,  100  U/ 
mL  penicillin/streptomycin,  and  4  mM  L-glutamine.  The  pSLIK-mtIDH  cell  lines 
(see  below)  were  maintained  as  above,  but  FBS  was  substituted  for  Tet-free 
FBS  (Clontech).  Cell  number  was  determined  using  an  automated  cell  counter 
(Nexcelom)  or  by  haemocytometer.  For  isotopic  labeling  experiments  in  the 
pSLIK-mtIDH  cell  lines,  cells  were  cultured  in  6-well  plates  in  glucose-  and 
glutamine-  free  DMEM,  supplemented  with  10%  dialyzed  Tet-free  FBS, 
100  U/mL  penicillin/streptomycin,  4  mM  L-glutamine,  and  10  or  15  mM  of 
the  appropriate  deuterated  glucose  tracer  ([3-^H,  95%  or  98%]glucose  or 
[4-^H,  94%  or  98%]glucose)  for  24  hr  or  48-72  hr  incubation,  respectively 
(Omicron  and  Cambridge  Isotope  Laboratories,  Inc.).  For  cholesterol  labeling 
experiments,  parental  HI  299  cells  were  cultured  in  DMEM  supplemented 
with  1%  FBS  for  two  passages  prior  to  72  hr  incubation  with  [3-^H]glucose. 
Cells  were  cultured  in  tracer  medium  for  24  hr  prior  to  the  addition  of  doxycy- 
cline  hyclate  (0.1  |ig/mL  in  water:  Sigma)  for24-48  hr  in  orderto  induce  mutant 
IDH  expression  and  accumulate  2HG.  Isotope-labeled  glycine  and  serine 
tracer  medium  was  prepared  from  custom  phenol  red-,  glucose-,  sodium 
pyruvate-,  amino  acid-,  and  sodium  bicarbonate-free  DMEM  (Hyclone  Labora¬ 
tories,  Inc.)  supplemented  with  10%  dialyzed  Tet-free  FBS,  3.7  g/L  sodium 
bicarbonate,  and  DMEM-levels  of  L-arginine,  L-cystine,  L-glutamine,  L-histi- 
dine,  L-isoleucine,  L-leucine,  L-lysine,  L-methionine,  L-phenylalanine,  L-thre- 
onine,  L-tryptophan,  L-tyrosine,  and  L-valine  prepared  as  a  lOOx  stock 
in  aqueous  acid  (pH  2.0).  A549  pSLIK  mtIDHI-C  and  mtlDH2-M  cells 
were  cultured  in  serine-  and  glycine-  free  DMEM  supplemented  with  either 
[2,3,3-^H3,  98%]serine  or  [3,3-^H2,  98%]serine  and  unlabeled  glycine 
(0.4  mM)  (Cambridge  Isotope  Laboratories,  Inc.),  or  [2,2-^H2,  98%]glycine 
and  unlabeled  serine  (0.4  mM)  (Cambridge  Isotope  Laboratories,  Inc.).  Cells 
were  cultured  in  the  presence  of  tracer  medium  for  24  hr  prior  to  doxycycline 
addition  (0.1  |.ig/mL)  for  a  further  48  hr. 

Generation  of  Cell  Lines  Stably  Expressing  Inducible  Forms  of 
FLAG-Tagged  Mutant  IDH 

To  generate  the  doxycycline-inducible  mutant  IDH  (mtIDH)  cell  lines,  full- 
length  cDNA  for  IDH1-R132H  and  IDH2-R172K  was  amplified  by  PCR  and 
cloned  into  the  p3xFLAG-CMV14  vector  (Sigma)  to  generate  C-terminal 
FLAG-tagged  constructs.  cDNA  for  IDH1-R132H-FLAG  and  IDH2-R172K- 
FLAG  was  then  amplified  by  PCR  and  cloned  into  the  pEN_TTmcs  entry  vector 


for  recombination  into  the  pSLIK-hygro  lentiviral  vector  (both  vectors  from 
Addgene  [Shin  et  al.,  2006]).  Lentiviruses  were  produced  by  transfecting 
HEK293T  cells  with  the  pSLIK-hygro-IDH1-R132H  or  pSLIK-hygro-IDH2- 
R172K  plasmids  along  with  the  lentiviral  packaging  plasmids  pMDLg/pRRE 
and  pRSV-Rev  and  the  envelope  plasmid  pMD2.G  (all  from  Addgene). 
Supernatants  containing  lentiviral  particles  were  collected  48  hr  after  transfec¬ 
tion  and  used  to  infect  subconfluent  H1299  and  A549  cells.  Infected  cells 
were  allowed  to  recover  for  24  hr  before  being  placed  under  selection  with 
350  i-ig/mL  hygromycin  (Invitrogen)  for  10  days.  Protein  expression  was 
induced  using  0.1  i.ig/mL  doxycycline  hyclate  (Sigma)  for  24-48  hr. 

Protein  Expression  Analysis  and  Cellular  Fractionation 

For  whole-cell  extracts,  cells  were  lysed  in  RIPA  buffer.  Mitochondrial  and 
cytoplasmic  fractions  were  prepared  as  previously  described  (Vender  Heiden 
et  al.,  1997).  Briefly,  cells  were  harvested  in  buffer  A  (250  mM  sucrose,  20  mM 
HEPES,  10  mM  KCI,  1.5  mM  MgCIs,  1  mM  EDTA,  1  mM  EGTA,  1  mM  DTT,  1  x 
EDTA-free  protease  inhibitor  cocktail  tablet  (Roche)  [pH  7.4])  and  broken  apart 
using  a  mechanical  homogenizer  (H  &  Y  Enterprise,  Redwood  City,  CA). 
Following  centrifugation  at  750  x  g  to  remove  unlysed  cells  and  nuclei,  mito¬ 
chondria  were  isolated  by  centrifuging  at  1 0,000  x  g  for  25  min.  The  resulting 
pellet  was  resuspended  in  buffer  A  and  represents  the  mitochondrial  fraction. 
The  remaining  supernatant  containing  cytoplasmic  and  membrane  proteins 
was  centrifuged  for  1  hr  at  100,000  x  g.  The  supernatant  from  this  final  spin 
represents  the  SI  00  fraction.  Protein  expression  was  analyzed  by  western 
blotting  using  antibodies  against  FLAG  (DYKDDDDK  tag,  Cell  Signaling), 
IDH1  (Santa  Cruz),  IDH2  (Abeam),  Cytochrome  C  (clone  7H8.2C12,  Abeam), 
Hsp  70  (Cell  Signaling). 

Metabolite  Extraction  and  GC-MS  Analysis 

Polar  metabolites  and  fatty  acids  were  extracted  using  methanol/water/chlo¬ 
roform  as  previously  described  (Metallo  et  al.,  2012).  Parental  and  pSLIK- 
mtlDH  cells  were  cultured  in  6-well  or  12-well  plates,  and  volumes  of  tracer 
media  and  extraction  buffers  were  adjusted  accordingly.  Derivatization  of 
both  polar  metabolites  and  fatty  acids  has  been  described  previously  (Metallo 
et  al.,  2012).  Briefly,  polar  metabolites  were  derivatized  to  form  methoxime- 
tBDMS  derivatives  by  incubation  with  2%  methoxylamine  hydrochloride  (MP 
Biomedicals)  in  pyridine  (or  MOX  reagent  (Thermo  Scientific)  followed  by  addi¬ 
tion  of  N-tert-butyIdimethylsilyl-N-methyltrifluoroacetamide  (MTBSTFA)  with 
1  %  tert-butyidimethylchlorosilane  (t-BDMCS)  (Regis  Technologies).  Nonpolar 
fractions,  including  triacylglycerides  and  phospholipids,  were  saponified  to 
free  fatty  acids  and  esterified  to  form  fatty  acid  methyl  esters  by  incubation 
with  2%  H2SO4  in  methanol.  Derivatized  samples  were  analyzed  by  GC-MS 
using  a  DB-35MS  column  (30  m  x  0.25mm  i.d.  x  0.25  |.im,  Agilent  J&W  Scien¬ 
tific)  installed  in  an  Agilent  7890A  gas  chromatograph  (GC)  interfaced  with  an 
Agilent  5975C  mass  spectrometer  (MS).  Mass  isotopomer  distributions  were 
determined  by  integrating  metabolite  ion  fragments  (Table  SI)  and  corrected 
for  natural  abundance  using  in-house  algorithms  adapted  from  Fernandez 
et  al.  (1996). 

Extraction  of  NAD+,  NADH,  NADP+,  and  NADPH  and  Analysis  by  LC- 
MS/MS 

The  extraction  protocol  for  NADPH  was  based  on  one  previously  described  by 
Fendt  et  al.  (201 3)  and  was  optimized  for  analysis  by  LC-MS/MS.  Briefly,  cells 
were  cultured  in  6-well  plates  over  the  course  of  30  min,  washed  once  in  ice- 
cold  water,  and  immediately  quenched  in  liquid  nitrogen.  A  total  of  200  |.il 
ice-cold  extraction  buffer  (40:40:20  acetonitrile/methanol/200  mM  NaCI, 
10  mM  Tris-HCI  [pH  9.2])  was  added  directly  to  the  cells.  Cells  were  scraped 
on  ice  and  cleared  by  centrifugation  at  4°C.  A  total  of  50  |.il  of  supernatant 
was  transferred  to  a  polypropylene  vial,  and  samples  were  analyzed  using  a 
Q  Exactive  Benchtop  LC-MS/MS  (Thermo  Fisher  Scientific). 

For  measurement  of  NAD(P)H  at  24  hr,  cells  were  cultured  in  10  cm  plates. 
After  24  hr  incubation  with  tracer,  approximately  1x10^  cells  were  washed  in 
ice-cold  0.9%  saline  and  immediately  quenched  in  1  ml  of  80%  methanol 
at  — 80°C.  Cells  were  scraped  on  dry  ice  and  cleared  by  centrifugation  at 
4°C.  Cleared  supernatant  was  transferred  to  Eppendorf  tube,  dried  under  vac¬ 
uum  using  a  CentriVap  (Labconco),  resuspended  in  water,  and  immediately 
loaded  onto  aXSELECTHSSXP  150  mm  x  2.1  mm  x  2.5  ^m  (Waters,  Milford, 
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MA)  with  an  UFLC  XR  HPLC  (Shimadzu,  Columbia,  MD)  coupled  to  an  AB 
SCIEX  Qtrap  5500  mass  spectrometer  (AB  SCIEX,  Framingham,  MA)  operating 
in  negative  ion  mode.  Mass  isotopomer  distributions  were  corrected  for  natu¬ 
ral  abundance  using  in-house  software  adapted  from  Fernandez  et  al.  (1996). 
Additional  information  regarding  chromatographic  separation,  mass  spec¬ 
trometry,  and  data  acquisition  can  be  found  in  the  Supplemental  Experimental 
Procedures. 

Isotopomer  Spectral  Analysis 

The  ISA  method  compares  a  measured  palmitate  mass  isotopomer  distribu¬ 
tion  to  one  that  is  simulated  using  a  reaction  network  for  palmitate  synthesis 
whereby  14  NADPH  molecules  are  consumed  to  form  one  palmitate  mole¬ 
cule.  Models  were  also  generated  for  myristate  and  stearate  synthesis 
whereby  12  or  16  NADPH  molecules  are  consumed  to  form  one  myristate 
or  stearate  molecule,  respectively.  Parameters  for  the  relative  enrichment  of 
the  lipogenic  NADPH  pool  from  a  given  pH]  tracer  and  the  percentage  of 
fatty  acids  that  are  de  novo  synthesized  are  extracted  from  a  best-fit  model 
using  the  INCA  MFA  software  package  (Figure  S2)  (Young,  2014).  The  95% 
confidence  intervals  for  both  parameters  were  estimated  by  evaluating  the 
sensitivity  of  the  sum  of  squared  residuals  between  measured  and  simulated 
palmitate  mass  isotopomer  distributions  to  small  flux  variations  (Antoniewicz 
et  al.,  2006). 

Ceil  Proliferation  Assays 

On  day  —1 , 1/10  of  a  confluent  10  cm  dish  of  cells  was  seeded  in  six  6-wells  of 
a  6-well  plate.  Twenty-four  hours  later,  cells  were  counted  on  an  automated 
cell  counter  (Nexcelom),  and  this  time  point  was  considered  Tq.  At  Tq,  all  other 
time  points  were  media  changed  to  ±  doxycycline  media  (0.1  iig/mL  doxycy- 
cline  hyclate  [Sigma]  in  water).  Cells  were  counted  every  24  hr  in  technical 
duplicate  and  biological  triplicate.  Media  was  changed  every  48  hr  to  prevent 
degradation  of  doxycycline  in  the  media. 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  five  figures,  one  table,  and  Supplemental 
Experimental  Procedures  and  can  be  found  with  this  article  at  http://dx.doi. 
org/10.1016/j.molcel.2014.05.008. 
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Abstract 

Oncogenic  mutations  in  isocitrate  dehydrogenase  1  and  2  (IDHl/2)  occur  in  several  types  of  cancer,  but  the 
metabolic  consequences  of  these  genetic  changes  are  not  fully  understood.  In  this  study,  we  performed  *®C 
metabolic  flux  analysis  on  a  panel  of  isogenic  cell  lines  containing  heterozygous  IDHl/2  mutations.  We  observed 
that  under  hypoxic  conditions,  IDHl-mutant  cells  exhibited  increased  oxidative  tricarboxylic  acid  metabolism 
along  with  decreased  reductive  glutamine  metabolism,  but  not  IDH2-mutant  cells.  However,  selective  inhibition 
of  mutant  IDHl  enzyme  function  could  not  reverse  the  defect  in  reductive  carboxylation  activity.  Furthermore, 
this  metabolic  reprogramming  increased  the  sensitivity  of  IDHl-mutant  cells  to  hypoxia  or  electron  transport 
chain  inhibition  in  vitro.  Lastly,  IDHl-mutant  cells  also  grew  poorly  as  subcutaneous  xenografts  within  a  hypoxic 
in  vivo  microenvironment.  Together,  our  results  suggest  therapeutic  opportunities  to  exploit  the  metabolic 
vulnerabilities  specific  to  IDHl  mutation.  Cancer  Res;  74(12);  1-15.  ©2014  AACR. 


Introduction 

Mutations  in  the  metabolic  enzymes  isocitrate  dehydroge¬ 
nase  1  and  2  (IDHl/2)  have  been  identified  in  a  variety  of  tumor 
types,  including  acute  myelogenous  leukemia  (AML),  gliomas, 
cholangiocarcinomas,  and  chondrosarcomas  (1-9).  These  muta¬ 
tions  are  almost  exclusively  heterozygous  point  mutations  that 
occur  in  specific  residues  within  the  cataljdic  pocket,  and  are 
suggestive  of  activating,  oncogenic  mutations.  Although  IDH 
mutants  are  no  longer  capable  of  efficiently  carrying  out  the 
normal  oxidative  reaction  [converting  isocitrate  and  NADP^  to 
a-ketoglutarate  (aKG),  CO2,  and  NADPH),  IDH  mutations  result 
in  a  novel  gain-of-function  involving  the  reductive,  NADPH- 
dependent  conversion  of  aKG  to  (Z))2-hydroxyglutarate  (2-HG; 
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refs.  10, 11).  2-HG  is  not  typically  present  at  high  levels  in  normal 
cells  but  accumulates  considerably  in  cells  with  IDHl/2  muta¬ 
tions  as  well  as  in  the  tumors  of  patients  with  IDHl/2  mutations, 
and  has  thus  been  termed  an  "oncometabolite  "  (10-12). 

Research  into  the  oncogenic  function  of  mutant  IDHl/2  has 
focused  in  large  part  on  the  effects  of  2-HG.  Numerous  reports 
have  linked  2-HG  accumulation  to  epigenetic  changes,  which 
are  thought  to  contribute  to  alterations  in  cellular  differenti¬ 
ation  status  (13-22).  Additional  mutant  IDH  phenotypes  have 
also  been  reported,  including  changes  in  collagen  maturation 
and  hypoxia  inducible  factor-la  (HIFla)  stabilization  (21, 23). 
These  changes  likely  occur  through  inhibition  of  aKG-depen- 
dent  dioxygenase  activity  by  high  levels  of  2-HG.  However,  the 
diverse  roles  that  aKG-dependent  dioxygenases  play  in  the  cell 
and  the  numerous  phenotypes  associated  with  mutant  IDH 
and  2-HG  suggest  that  the  phenotypes  downstream  of  2-HG 
induction  could  be  cell  type-  or  context-specific.  We  hypoth¬ 
esize  that  metabolic  alterations  induced  by  IDH  mutations 
may  also  be  present  and  might  be  a  general  phenotype  that 
offers  additional  approaches  to  target  these  tumors.  Previous 
work  suggests  that  overexpression  of  mutant  IDH  alters  the 
levels  of  several  metabolites  (24)  and  leads  to  increased  sen¬ 
sitivity  to  glutaminase  inhibitors  (25).  Studies  by  Leonardi  and 
colleagues  have  indicated  that  the  IDHl-mutant  enzymes 
compromise  the  ability  of  this  enzyme  to  catalyze  the  reductive 
carboxylation  reaction  (26).  However,  it  is  unclear  how 
IDH  mutations  affect  central  carbon  metabolism  in  the  het¬ 
erozygous  cellular  setting,  and  further  exploration  into  how 
these  metabolic  differences  could  be  therapeutically  exploited 
is  warranted.  An  important  distinction  between  IDHl  and 
IDH2  is  their  localization  in  the  cytosol/peroxisome  and 
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mitochondria,  respectively,  which  may  influence  the  ultimate 
metabolic  phenotype  of  tumor  cells  with  mutations  in  either 
enzyme. 

Systems-based  approaches  that  use  stable  isotope  tracers  (e.g., 
[^^CJglutamine),  mass  spectrometry,  and  network  modeling 
to  estimate  metabolic  fluxes  offer  a  unique  means  of  character¬ 
izing  intracellular  metabolism  (27).  To  understand  the  metabolic 
impact  of  heterozygous  IDH  mutation  in  vitro,  we  have  appUed 
metabohc  flux  analysis  (MFA)  to  a  panel  of  cell  lines  that 
differ  only  with  respect  to  their  IDHl-  and  lDH2-mutant  status. 
Using  this  approach,  we  have  characterized  how  cells  with  wild- 
type  (WT)  and  mutant  lDHl/2  respond  to  hypoxia  and  phar¬ 
macologic  induction  of  mitochondrial  dysfunction. 

Materials  and  Methods 

Cell  culture 

HCT116  and  MCF-lOA  isogenic  clones  were  obtained  from 
Horizon  Discovery  Ltd  and  IDH  1/2  mutational  status  was 
verified  by  sequencing  (28).  HCT116  cells  were  cultured  in 
McCoy's  5A  Modified  Medium  with  10%  FBS.  D-2-HG  treat¬ 
ments  were  done  at  10  mmol/L  and  replenished  every  48  hours. 
MCF-lOA  cells  were  cultured  as  described  previously  (28).  HT- 
1080,  SW1353,  A549,  and  143B  cells  were  cultured  in  DMEM 
supplemented  with  10%  FBS.  HT-1080  and  SW1353  cells  were 
obtained  from  the  ATCC,  and  cells  were  tested  and  authenti¬ 
cated  by  single-nucleotide  polymorphism  fingerprinting.  A549 
cells  were  obtained  from  ATCC  and  were  not  further  tested  or 
authenticated.  143B  cells  were  kindly  provided  by  Dr.  Leonard 
Guarente  (Massachusetts  Institute  of  Technology,  Cambridge, 
MA)  and  were  not  further  tested  or  authenticated.  Cells  were 
routinely  cultured  in  normoxia  (21%  O2)  and  then  moved  to 
hypoxia  (l%-3%  O2,  as  indicated  in  the  figure  legends)  for 
48  to  72  hours  where  indicated.  Generation  of  the  p°  cells  is 
described  in  the  Supplementary  Methods.  Xenograft  assays  are 
described  in  Supplementary  Methods. 

Steady-state  labeling  of  organic,  amino,  and  fatty  acids  was 
accomplished  by  culturing  subconfluent  cells  in  triplicate  in 
tracer  medium  for  72  hours  in  a  6-well  plate.  Labeling  studies  of 
HCT 116,  SW1353,  HT 1080,  A549,  and  143B  cells  were  performed 
in  glucose  or  glutamine-free  DMEM  containing  10%  FBS  and 
17.5  mmol/L  [l,2-^^C2]glucose,  4  mmol/L  [U-’^CsJglutamine,  or 
4  mmol/L  [5-*^C]glutamine.  For  the  HCT116  isogenic  cells,  the 
initial  seeding  density  was  150,000  cells  per  well  except  for  the 
IDHl  R132H/-f  2H1,  IDHl  R132C/+  2A9,  and  IDHl  R32C/+ 
3A4,  which  were  250,000  cells  per  well.  Labeling  studies  of  the 
MCF-lOA  cells  were  done  in  glutamine-free  DMEM  containing 
4  mmol.mg/L  [U-^^Cslglutamine,  5%  horse  serum,  20  ng/mL 


EGF,  10  pg/mL  insuhn,  0.5  pg/mL  hydrocortisone,  and 
100  ng/ mL  cholera  toxin.  Gas  chromatography  mass  spectrom¬ 
etry  (GC-MS)  analysis  is  described  in  Supplementary  Methods. 

MFA 

^^C  MFA  was  conducted  using  INCA,  a  software  package 
based  on  the  EMU  framework  (http://mfa.vueinnovations. 
com;  ref  29).  Intracellular  concentrations  of  free  metabolites 
and  intra-  and  extracellular  fluxes  were  assumed  to  be  constant 
over  the  course  of  the  tracing  experiment.  Fluxes  through  a 
metabolic  network  comprising  of  glycolysis,  the  pentose  phos¬ 
phate  pathway,  the  TCA  cycle,  biomass  synthesis,  and  fatty  acid 
synthesis  were  estimated  by  minimizing  the  sum  of  squared 
residuals  between  experimental  and  simulated  MIDs  and 
extracellular  fluxes  using  nonlinear  least  squares  regression 
(30).  The  best  global  fit  was  found  after  estimating  100  times 
using  random  initial  guesses  for  all  fluxes  in  the  network.  A 
statistical  test  was  applied  to  assess  the  goodness-of-fit  using  a 
of  0.01.  The  95%  confidence  intervals  for  all  fluxes  in  the 
network  were  estimated  by  evaluating  the  sensitivity  of  the 
sum  of  squared  residuals  to  flux  variations  (30).  Isotopomer 
Spectral  Analysis  was  performed  as  previously  described  (31). 
See  Supplementary  Methods  for  further  details  of  MFA. 

Reagents 

The  following  reagents  were  used  at  the  doses  indicated  and 
as  described  in  the  text/figure  legends:  [l,2-^^C2]glucose, 
[3-^^C]glucose,  [U-^^Cs]glutamine,  [l-^^C]glutamine,  and 
[5-^^C]glutamine  (all  from  Cambridge  Isotope  Laboratories); 
IDH-C277  (Xcessbio);  HIFla  antibody  (610958,  BD  Bioscinces). 
Synthesis  of  IDHli  A  is  described  in  Supplementary  Methods. 

Determination  of  oxygen  consumption 

HCT  116  cells  were  grown  at  either  normoxia  or  hypoxia  (3% 
O2),  and  respiration  was  measured  using  an  XF®96  analyzer 
(Seahorse  Bioscience).  Cell  growth  and  assays  at  3%  O2  were 
conducted  using  the  Coy  Dual  Hypoxic  Chambers  for  Seahorse 
XF®  Analyzer  (Coy  Laboratory  Products,  Inc.)  as  described  in 
Supplementary  Methods. 

Proliferation  assays 

To  calculate  doubling  time,  cells  were  trypsinized  and  viable 
cells  were  quantified  on  a  ViCell  (Beckman-Coulter).  Doubling 
times  are  presented  as  the  average  of  three  or  more  indepen¬ 
dent  experiments. 

To  generate  longer-term  growth  curves,  cells  were  plated  at 
3,000  cells  per  well  in  a  96-well  plate  in  triplicate.  Twenty-four 


Figure  1 .  Isogenic  IDH1  mutation  compromises  metabolic  reprogramming  under  hypoxia.  A,  carbon  atom  (represented  by  circles)  transitions  and  tracers  used 
to  detect  changes  influx:  [U-^^Cslglutamine  (purple  circles)  or  [5-^^C]glutamine  (circle  with  red).  The  fifth  carbon  is  lost  during  oxidative  TCA  metabolism  but  is 
retained  on  citrate,  AcCoA,  and  palmitate  in  the  reductive  pathway  (green  arrows).  B,  citrate  MID  labeling  from  [U-^^Cslglutamine  from  HCT116  parental 
and  HCT1 16  IDH1  R132H/+  clone  2H1  cells  cultured  in  normoxia  or  hypoxia  (2%  oxygen)  for  72  hours.  Data,  representative  of  more  than  three  independent 
experiments.  Inset  highlights  changes  in  %  M5  citrate.  C,  contribution  of  [5-^^C]glutamine  to  lipogenic  AcCoA  from  cells  cultured  as  in  Fig.  IB.  D, 
uptake  and  secretion  fluxes  for  cells  cultured  as  in  Fig.  1 B.  E-H,  a-Ketoglutarate  dehydrogenase  (E),  reductive  IDH  (E),  oxidative  IDH  (G),  and  glutamine 
anaplerosis  (H)  flux  estimates  and  95%  confidence  intervals  by  the  MFA  model.  I,  graphical  representation  of  fluxes  determined  via  MFA.  Arrow  thickness, 

level  of  flux  in  HCT  116  cells  cultured  in  hypoxia.  Color,  fold  difference  between  hypoxic  and  normoxic  parental  HCT  116  cells  (left)  or  between  hypoxic  HCT  116 
IDHl  R132H/+  2H1  cells  and  hypoxic  HCT1 16  parental  cells  (right).  *,  metabolites  that  were  modeled  as  existing  in  separate  mitochondrial  and  cytosolic 
pools.  aKG,  a-Ketoglutarate;  cit,  citrate;  fum,  fumarate;  glue,  glucose;  glu,  glutamate;  gin,  glutamine;  lac,  lactate;  mal,  malate;  oac,  oxaloacetate; 
pyr,  pyruvate.  See  also  Supplementary  Methods,  Supplementary  Fig.  S2,  and  Supplementary  Tables  S1-S4  for  details  of  MFA  model,  results,  and  data. 
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Figure  2.  Compromised  reductive  TCA  metabolism  is  specific  to  cells  with  mutant  IDH1 .  A,  relative  level  of  reductive  glutamine  metabolism,  determined  by  M5 
labeling  ofcitratefrom[U-^^C5]glutamineinHCT11 6  parental  cells  with  or  with  out  lOmmol/L  D-2-HG,orHCT116IDH1/2-mutant  isogenic  cells  cultured  in  normoxia 
or  hypoxia  (2%  oxygen).  Percentage  of  M5  citrate  levels  are  normalized  to  HCT 116  parental  cells  in  normoxiafor  each  experiment.  {Continued  on  the  following  page.) 
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hours  later,  the  indicated  treatment  was  started  and  con- 
fluency  measurements  were  taken  every  12  hours  for  108  to 
216  hours  using  an  Incucjde  Kinetic  Imaging  System  (Essen 
BioScience).  Confluency  data  were  modeled  using  a  general¬ 
ized  logistic  growth  equation  (equation  1),  and  the  maximum 
growth  rate  was  estimated  using  nonlinear  regression: 


where  U  and  L  represent  upper  and  lower  asymptotes,  to 
represents  the  time  at  which  cell  confluency  reaches  50%,  and 
Mmax  represents  the  maximum  growth  rate  per  hour. 

Pharmacologic  profiling  of  the  CCLE  was  performed  as 
previously  described  (32).  The  growth  inhibition  assays  are 
described  in  Supplementary  Methods. 

Statistical  analysis 

AU  results  shown  as  averages  of  multiple  independent 
experiments  are  presented  as  mean  ±  SE;  results  shown  as 
averages  of  technical  replicates  are  presented  as  mean  ±  SD. 
P  values  were  calculated  using  a  Student  two-tailed  t  test;  *, 
P  value  between  0.005  and  0.05;  **,P  value  between  0.001  and 
0.005;  ***,  P  value  <0.001.  All  errors  associated  with  MFA  and 
ISA  of  lipogenesis  are  95%  confidence  intervals  determined  via 
sensitivity  analysis. 

Results 

Mutant  IDHl  compromises  metabolic  reprogramming 
under  hypoxia 

We  and  others  have  previously  demonstrated  that  tricar¬ 
boxylic  acid  (TCA)  metabolism  is  reprogrammed  under  hyp¬ 
oxia,  and  flux  through  WT  IDHl  and/or  IDH2  become  critical 
in  these  contexts  (31,  33-37).  Oncogenic  mutations  in  IDHl 
and  IDH2  mitigate  these  enzymes'  WT  function  and,  in  par¬ 
ticular,  the  ability  to  catalyze  the  reductive  carboxylation 
reaction  while  inducing  a  neomorphic  activity  that  results  in 
the  accumulation  of  D-2-HG  (11,  26).  Therefore,  we  hypothe¬ 
sized  that  cancer  cells  harboring  mutations  in  either  IDHl  or 
IDH2  may  be  compromised  in  their  ability  to  modulate  metab- 
ohsm  under  low  oxygen  tensions.  To  identify  metabohc  liabil¬ 
ities  induced  by  IDHl  mutations,  we  applied  MEA  to 
isogenic  HCTl  16  colon  cancer  cells  with  WT  IDHl/2  (parental) 
or  a  heterozygous  IDHl  mutation,  IDHl  R132H/ -f  (clone  2H1). 
The  IDHl-mutant,  but  not  WT,  cells  have  previously  been 
shown  to  produce  high  levels  of  2-HG  (28).  To  gauge  relative 
flux  through  the  TCA  cycle,  each  cell  line  was  cultured  in  the 
presence  of  [U-^^Cs]glutamine  (uniformly  ^^C-labeled  gluta¬ 
mine)  under  normoxic  or  hypoxic  (2%  oxygen)  conditions  for 


72  hours,  and  isotope  enrichment  in  various  metabolites  was 
determined  via  mass  spectrometry  (Fig.  lA).  Both  cell  lines 
displayed  decreased  oxidative  TCA  flux  (as  evidenced  by 
decreased  M3  aKG)  in  hypoxia  (Supplementary  Fig.  SIA). 
Although  minimal  changes  in  labeling  were  detected  when 
comparing  the  mass  isotopomer  distribution  (MID)  of  citrate 
in  each  cell  type  grown  in  normoxia  (designated  as  MO,  Ml,  M2, 
etc.  mass  isotopomers,  corresponding  to  ion  fragments  con¬ 
taining  zero,  one,  or  two  ^^C,  respectively),  more  significant 
deviations  occurred  in  cells  proliferating  under  hypoxia  (Fig. 
IB).  Parental  cells  under  hypoxia  exhibited  increased  M5 
labeling  indicative  of  reductive  carboxylation  (Fig.  IB),  as  has 
been  seen  previously  for  many  WT  IDHl/2  cell  lines  (31, 33, 35). 
In  contrast,  the  IDHl  R132H/-f  cells  showed  only  a  slight 
increase  in  the  abundance  of  this  mass  isotopomer  under 
hypoxia  (Fig.  IB).  M5  citrate  can  also  be  produced  via  M6 
citrate  and  glutaminolysis  (38);  however,  no  increase  in  the  low 
basal  levels  of  M6  citrate  was  observed  in  hypoxia  (Fig.  IB).  We 
observed  similar  changes  in  the  labeling  of  other  TCA  meta¬ 
bolites,  including  M3  fumarate,  malate,  and  aspartate  (derived 
from  oxaloacetate),  further  supporting  our  finding  that  IDHl- 
mutant  cells  display  compromised  reductive  glutamine  metab¬ 
olism  in  hypoxia  (Fig.  lA  and  Supplementary  Fig.  SlB-SlD).To 
further  quantify  the  metabolism  of  glutamine  through  the 
reductive  carboxylation  pathway  in  these  cells,  we  determined 
the  contribution  of  [5-'^C]glutamine  to  palmitate  synthesis 
using  isotopomer  spectral  analysis  (ISA),  as  this  tracer  specif¬ 
ically  labels  acetyl  coenzyme  A  (AcCoA)  through  the  reductive 
carboxylation  pathway  (Fig.  lA;  ref  39).  Although  parental  cells 
increased  the  contribution  of  glutamine  to  lipogenic  AcCoA 
almost  5-fold  under  hypoxia,  cells  with  a  mutant  IDHl  allele 
were  compromised  in  their  ability  to  increase  this  reductive 
carboxylation  flux  (Fig.  1C). 

To  characterize  the  metabolic  phenotype  of  HCT116  cells 
with  WT  or  IDHl  R132H/-f  in  a  more  comprehensive  and 
unbiased  manner,  we  incorporated  uptake/secretion  fluxes 
and  mass  isotopomer  data  into  a  network  of  central  carbon 
metabolism.  This  model  included  glycolysis,  the  pentose  phos¬ 
phate  pathway  (PPP),  TCA  metabolism,  and  various  biosyn¬ 
thetic  fluxes  using  [U-^^Cslglutamine  and  [1,2- ^^C2] glucose  (for 
the  oxidative  PPP  bifmcation),  as  these  tracers  provide  optimal 
flux  resolution  throughout  central  carbon  metabolism  (40).  An 
elementary  metabolite  unit  (EMU)-based  algorithm  was  used 
to  estimate  fluxes  and  associated  confidence  intervals  in  the 
network  (41,  30),  and  a  detailed  description  of  the  model, 
assumptions,  and  the  complete  data  set  are  included  as 
Supplementary  Material.  As  expected,  glucose  and  lactate 
fluxes  were  significantly  increased  by  hypoxia  in  both  cell 
lines,  and  significant  2-HG  secretion  occurred  only  in  IDHl 


{Continued.)  B,  Western  blot  of  HCT116  parental,  HCT116  IDHl  R132H/+  clone  2C11,  and  HCT116  IDHl  R132H/+  clone  2H1  showing  levels  of  IDH1 
R132H  and  total  levels  of  IDHl.  C,  ratio  of  a-ketoglutarate  to  citrate  from  cells  cultured  as  In  Fig.  2A.  D,  contribution  of  [U-^^Cslglutamlne  to  lipogenic 
AcCoA  from  cells  cultured  as  in  Fig.  2A.  Four  independent  experiments  are  shown.  E,  correlation  of  reductive  glutamine  metabolism  in  2%  oxygen 
(as  measured  by  %  M5  citrate  from  [U-'^Cslglutamine)  to  2-HG  ng/1  million  cells.  Note  that  HCT1 16  IDH2  R172K/+  45G9  and  59F10  are  not  included  in 
this  figure.  F,  relative  level  of  reductive  glutamine  metabolism,  determined  by  M5  labeling  of  citrate  from  [U-^^Cslglutamine  in  HCTl  1 6  parental  cells,  which 
were  stably  infected  with  doxycycline-inducible  shRNA's  targeting  IDHl  or  IDH2,  or  nontargeting  control  (NTC).  Cells  were  cultured  for  48  hours  in 
normoxia  or  hypoxia  (1%  oxygen)  in  the  presence  of  100  ng/mL  doxycycline.  Citrate  labeling  was  determined  via  liquid  chromatography/mass 
spectrometry.  G,  relative  level  of  reductive  glutamine  metabolism,  determined  by  M5  labeling  of  citrate  from  [U-'^Cslglutamine  in  MCF-I OA  cells  cultured 
in  normoxia  or  hypoxia  (2%  oxygen)  for  72  hours.  H,  contribution  of  [U-’^Cslglutamine  to  lipogenic  AcCoA  from  cells  cultured  as  in  Fig.  2G. 
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R132H/+  cells  (Fig.  ID).  Notably,  2-HG  secretion  was  elevated 
under  hypoxia,  consistent  with  previous  observations  of  2-HG 
accumulation  at  low  oxygen  tension  (33). 

The  modeling  data  comparing  parental  HCT116  cells  grown 
in  normoxia  and  hypoxia  highlight  some  of  the  important 
metabolic  changes  that  occur  at  low  oxygen  tensions  (Fig.  lE-I, 
Supplementary  Fig.  S2,  and  Supplementary  Tables  S1-S4).  In 
the  HCT116  parental  cells,  pyruvate  dehydrogenase  and  oxi¬ 
dative  TCA  metabolism  are  decreased  under  hypoxia,  whereas 
pyruvate  cycling  through  malic  enzyme  (ME)  and  pyruvate 
carboxylase  (PC)  are  elevated  under  these  conditions.  Parental 
cells  increased  reductive  IDH  flux  several  fold,  such  that  net 
IDH  flux  slightly  favored  the  direction  of  reductive  carboxyl- 
ation  (Fig.  IF,  G,  and  I).  As  with  the  MID  changes  (Fig.  IB),  only 
modest  changes  in  metabolism  were  detected  when  compar¬ 
ing  parental  HCT116  cells  with  IDHl  R132H/-f  2H1  grown 
under  normoxia  (Fig.  lE-H  and  Supplementary  Tables  S1-S4). 
However,  mutant  cells  maintained  high  oxidative  IDH  and 
aKG-dehydrogenase  (aKGDH)  fluxes  and  were  unable  to 
induce  reductive  carboxylation  under  hypoxia  relative  to  the 
parental  cells  (Fig.  lE-G,  I,  and  Supplementary  Tables  S1-S4). 
This  oxidative  TCA  flux  was  maintained  by  increased  gluta¬ 
mine  anaplerosis  and  flux  through  ME  and  PC  (Fig.  IH  and  I 
and  Supplementary  Tables  S1-S4).  Overall,  these  results  dem¬ 
onstrate  that  significant  reprogramming  of  TCA  metabolism 
occurs  in  cells  at  2%  oxygen,  and  expression  of  IDHl  R132H/-f 
abrogates  the  cells'  ability  to  respond  appropriately  to  hypoxic 
stress. 

Compromised  reductive  TCA  metabolism  is  specific  to 
cells  with  mutant  IDHl 

The  MFA  results  above  suggest  that  heterozygous  IDHl 
mutations  disrupt  the  metabolic  response  of  cells  to  hypoxic 
stress.  To  determine  whether  this  metabolic  deficiency  is 
common  to  cells  with  either  IDHl  or  IDH2  mutations,  we 
interrogated  a  panel  of  previously  reported  IDHl-  and  IDH2- 
mutant  isogenic  cells  (28)  and  measured  the  ability  of  each  to 
initiate  reductive  carboxylation  under  2%  oxygen.  With  the 
exception  of  IDH2  R140Q/-f  cells,  these  mutant  cell  lines  exhibit 
a  >25-fold  increase  in  2-HG  compared  with  parental  HCT116 
cells  (28).  We  cultured  each  clone  with  [U-^^CsJglutamine  under 
normoxia  and  hypoxia,  quantifying  M5  citrate  abundance 
(see  Fig.  lA)  to  determine  the  relative  extent  of  reductive 
carboxylation  induction.  All  but  one  of  the  IDHl-mutant  clones 
were  consistently  compromised  in  their  ability  to  increase 
reductive  carboxylation  activity  under  hypoxia  (Fig.  2A).  The 
one  exception  being  the  IDHl  R132H/+  2C11  clone,  which 
showed  a  weaker  phenotype  relative  to  the  other  IDHl-mutant 
clones  (Fig.  2A).  This  is  likely  explained  by  a  lower  level  of  IDHl 
R132H  protein  (Fig.  2B)  and  lower  level  of  2-HG  than  the  IDHl 
R132H/-f  2H1  clone  (28). 

Unlike  mutant  IDHl  cells,  HCT116  cells  with  IDH2  muta¬ 
tions  exhibited  levels  of  M5  citrate  under  hypoxia  that  were 
comparable  with  the  parental  cells  (Fig.  2A).  Addition  of  10 
mmol/L  D-2-HG  to  the  culture  media  also  had  a  minimal  effect 
on  TCA  metabolism  (Fig.  2A).  These  data  suggest  that  high  2- 
HG  levels  alone  are  unable  to  inhibit  reductive  carboxylation 
activity,  even  though  this  dose  of  exogenous  D-2-HG  is  suffi- 


Figure  3.  Cells  with  endogenous  IDH1  and  IDH2  mutations  respond 
differently  to  mitochondrial  stress.  A,  relative  level  of  reductive  glutamine 
metabolism,  determined  by  M5  labeling  of  citrate  from  [U-^^Cslglutamine 
in  HT-1080  and  SW1353  cells  cultured  in  normoxia  or  hypoxia  (2% 
oxygen)  for  72  hours.  B,  contribution  of  [5-^^C]glutamine  to  lipogenic 
AcCoA  from  HT-1 080  and  SW1 353  cells  cultured  in  normoxia  or  hypoxia 
(1  %  oxygen)  for  72  hours.  C,  relative  level  of  reductive  glutamine 
metabolism,  determined  by  M5  labeling  of  citrate  from  [U-^^Cslglutamine 
in  HT-1080  and  SW1353  WT  or  p°  cells.  D,  contribution  of 
[U-^^Cslglutamine  and  [U-'^Celglucose  to  lipogenic  AcCoA  from  cells 
cultured  as  in  Fig.  3C.  E,  relative  level  of  reductive  glutamine  metabolism, 
determined  by  M5  labeling  of  citrate  from  [U-^^Cslglutamine  in  A549  and 
143B  WT  or  p°  cells.  F,  contribution  of  [U-^^Cslglutamine  to  lipogenic 
AcCoA  from  cells  cultured  as  in  Fig.  3E. 
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Figure  4.  Mutant  IDH1  affects  TCA 
metabolism  in  vivo.  A  and  B, 
box-whisker  plots  showing 
[1-^^C]glutamine  enrichment  in 
plasma  (A)  and  tumor  (B)  from  mice 
with  xenografts  derived  from  the 
indicated  cell  lines.  Black  dot, 
mean;  center  line,  median;  box, 
interquartile  range;  whiskers, 
maximum  and  minimum  values. 

C,  relative  2-HG  levels  from 
xenografts,  normalized  to  HCT116 
parental  xenografts.  D,  percentage 
of  Ml  aKG  derived  from 
[1-^^C]glutamine  in  tumor 
xenografts.  E,  nMols  aKG/nMols 
citrate  ratio  from  xenografts  grown 
from  the  indicated  HCT116  cells. 

F,  glutamine  anaplerosis,  as 
determined  by  M1  labeling  of  aKG 
relative  to  M1  labeling  of  glutamine 
from  xenografts  grown  from  the 
indicated  HCT116  cells. 
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cient  to  induce  the  2-HG-dependent  epithelial-mesenchymal 
transition  (EMT)  phenotype  previously  described  in  these 
cell  lines  (28).  Consistent  changes  were  observed  when 
normalizing  M5  citrate  to  M5  glutamate  (Supplementary 
Fig.  S3A),  and  levels  of  M6  citrate  were  not  high  in  any  of  the 
clones  (Supplementary  Fig.  S3B),  providing  evidence  that 
such  changes  are  specific  to  the  IDH/aconitase  node  of 
metabolism.  Similar  trends  were  also  observed  when  mea¬ 
suring  the  ratio  of  aKG/ citrate  under  normoxia  and  hypoxia 
(Fig.  2C),  another  metric  that  describes  the  extent  of  reduc¬ 
tive  versus  oxidative  IDH  flux  (36,  42).  Finally,  the  contri¬ 
bution  of  glutamine  to  lipogenic  AcCoA  under  hypoxia  was 
significantly  lower  in  cells  with  IDHl  mutations  but  not 
those  with  IDH2  mutations  or  exogenous  2-HG  (Fig.  2D). 
Overall,  the  extent  that  each  IDHf-mutant  cell  line  produced 
2-HG  correlated  with  their  ability  to  activate  reductive 
carboxylation  flux  under  hypoxia  (Fig.  2E),  whereas  IDH2- 
mutant  cells  did  not  adhere  to  this  trend. 

We  next  conducted  shRNA-mediated  knockdown  of  IDHl 
and  IDH2  in  parental  HCT116  cells  to  examine  the  roles  of  WT 


IDHl  and  IDH2  in  mediating  reductive  glutamine  metabolism 
(Supplementary  Fig.  S3C  and  S3D).  Knockdown  of  IDHl 
decreased  levels  of  M5  citrate  in  cell  populations  cultured 
with  [U-^^CsJglutamine,  whereas  cells  expressing  shRNAs  tar¬ 
geting  IDH2  exhibited  the  same  or  higher  M5  citrate  (Fig.  2F). 
These  results  are  consistent  with  previous  studies  in  other  cell 
lines  (31),  highhghting  the  importance  of  WT  IDHl  expression 
in  reprogramming  TCA  metabolism  under  hypoxia,  and  fur¬ 
ther  suggest  that  IDHl  mutation  selectively  impedes  WT  IDHl 
function  in  these  cells. 

To  determine  whether  this  mutant  IDHl-induced  metabolic 
deficiency  is  specific  to  the  HCT116  genetic  background  or 
more  broadly  applicable  to  cells  of  different  tissue  origin,  we 
performed  similar  analyses  using  MCF-lOA  immortalized 
mammary  epithehal  cells  with  heterozygous  IDHl  mutations 
(28).  When  cultured  for  3  days  under  normoxia  or  hypoxia,  two 
distinct  IDHl  R132H/-f  clones  were  compromised  in  their 
ability  to  generate  M5  citrate  or  lipogenic  AcCoA  from 
[U-^^Cslglutamine  (Fig.  2G  and  H  and  Supplementary  Fig. 
S3E).  Thus,  cells  with  heterozygous  IDHl  mutations,  but  not 
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Figure  5.  Inhibition  of  mutant  IDH1  does  not  rescue  reprogramming  of  TCA  metabolism.  A,  structure  and  biochemical  data  for  mutant  IDH1  inhibitor, 
IDH1  i  A.  Note  that  this  is  the  (S)  enantiomer.  B,  2-HG  levels  in  HCT1 1 6  IDH1  R1 32H/+  2H1  and  HCT1 1 6  IDH1  R1 32C/+  2A9  cell  lines  treated  with  the 
indicated  concentrations  of  IDH1  i  A  for  3  days.  C,  doubling  time  of  cells  cultured  with  or  without  10  |imol/L  of  IDHIi  A.  D,  change  in  total  metabolite  levels  of 
HCT116  IDH1  R132H/+ 2H1  and  HCT116IDH1  R132C/+ 2A9  cells  cultured  in  the  presence  or  absence  of  10  pmol/L  of  IDH1  i  A  for  3  or  31  days,  the 
final  72  hours  of  which  the  cells  are  grown  in  hypoxia  (2%  oxygen).  E,  citrate  MID  labeling  from  [U-^^CsJglutaminefrom  cells  cultured  as  in  Fig.  5D.  Dotted  line 
over  M5  citrate  represents  average  %  M5  citrate  observed  in  HCT1 16  parental  cells  cultured  in  hypoxia.  {Continued  on  the  following  page.) 
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IDH2  mutations  or  exogenous  2-HG,  are  compromised  for 
glutamine  reductive  carboxylation  under  hypoxia. 

Cells  with  endogenous  IDHl  and  IDH2  mutations 
respond  differently  to  mitochondrial  stress 

To  examine  whether  these  trends  are  observed  in  cell 
lines  from  cancers  with  endogenous  IDH  mutations,  we 
evaluated  the  ability  of  two  cell  lines  harboring  IDHl  or 
IDH2  mutations  to  activate  reductive  carboxylation  under 
conditions  of  mitochondrial  stress.  When  switched  to  hyp¬ 
oxic  growth,  HT1080  IDHl  R132C/+  fibrosarcoma  cells 
exhibited  a  significantly  decreased  ability  to  induce  reduc¬ 
tive  glutamine  metabolism  in  comparison  with  SW1353 
IDH2  R172S/-f  chondrosarcoma  cells  (Fig.  3A  and  Supple¬ 
mentary  Fig.  S4A).  HT-1080  cells  also  used  less  glutamine  for 
de  novo  lipogenesis  than  the  SW1353  cells  in  1%  oxygen 
tension  (Fig.  3B).  Thus,  a  cell  line  with  an  endogenous 
mutation  in  IDHl,  hut  not  IDH2,  displays  compromised 
reductive  glutamine  metabolism  in  hypoxia. 

In  addition  to  low  oxygen  tensions,  an  alternative  means 
of  inducing  reductive  TCA  metabolism  is  through  the  inhi¬ 
bition  of  oxidative  phosphorylation  (OXPHOS;  refs.  34,  42). 
To  compromise  OXPHOS,  we  generated  p“  cells  that  lack  a 
functional  electron  transport  chain  (ETC)  from  various  cell 
lines  using  established  methods  (43).  As  expected,  oxidative 
mitochondrial  metabolism  was  virtually  extinguished,  as 
evidenced  by  M3  and  M3/M5  labeling  of  aKG  in  p“  cells 
generated  from  IDH-mutant  cells  (HT1080  and  SW1353)  or 
other  cancer  cell  lines  with  WT  IDHl/2  (143B  osteosarcoma, 
A549  non-small  cell  lung  cancer;  Supplementary  Fig.  S4B- 
S4E).  However,  HT1080  IDHl  R132C/-f  p”  cells  were  com¬ 
promised  in  their  ability  to  generate  citrate  and  palmitate 
through  reductive  glutamine  metabolism,  whereas  SW1353 
and  both  of  the  IDHl/2  WT  p“  cell  lines  were  able  to 
efficiently  induce  reductive  carboxylation  and  use  glutamine 
for  lipid  synthesis  (Fig.  3C-F).  Both  the  HT1080  and  SW1353 
p“  cells  continued  to  use  glucose  for  lipid  synthesis,  and  this 
contribution  was  higher  in  the  HT1080  p**  cells  (Fig.  3D).  The 
increased  glucose  utilization  in  HT1080  p“  cells  compared 
with  SW1353  p“  cells  was  facilitated  by  anaplerosis  through 
pyruvate  carboxylase,  as  demonstrated  by  increased  labeling 
in  TCA  intermediates  from  [3-^^C]glucose  (Supplementary 
Fig.  S4F  and  S4G).  These  results  provide  evidence  that 
hypoxia  and  mitochondrial  dysfunction  drive  reprogram¬ 
ming  of  the  TCA  cycle,  and  cells  with  spontaneously 
acquired  IDHl  mutations  are  unable  to  efficiently  reprogram 
metabolism  to  induce  reductive  glutamine  carboxylation. 

Mutant  IDHl  affects  TCA  metabolism  in  vivo 

The  metabolic  deficiencies  of  IDHl-mutant  cells  occur  at 
oxygen  tensions  that  are  likely  to  occur  in  solid  tumors  and 
some  normal  tissues  (44).  To  determine  whether  these 


metabolic  phenotypes  arise  in  vivo,  we  generated  subcuta¬ 
neous  xenografts  using  parental,  IDHl  R132H/-f  2H1,  and 
IDHl  R132C/-f  2A9  HCT116  cells.  After  tumors  achieved  a 
minimum  diameter  of  0.8  cm,  mice  were  infused  with 
[l-^^C]glutamine  for  6  hours  to  achieve  steady-state  isotope 
enrichment  in  plasma  and  tumor  (Fig.  4A  and  B;  ref.  45).  A 
targeted  metabolomic  analysis  was  performed  on  plasma 
and  tumor  extracts  to  quantify  isotope  enrichment  and 
metabolite  abundances.  Significant  2-HG  was  detected  only 
in  the  IDHl-mutant  tumors  (Fig.  4C).  Insufficient  isotope 
enrichment  was  achieved  in  plasma  and  intratumoral  glu- 
tamine/aKG  to  detect  label  on  citrate  via  reductive  carbox¬ 
ylation  (Fig.  4A,  B,  and  D).  However,  in  agreement  with  the 
results  obtained  from  in  vitro  studies,  the  aKG/citrate  ratio 
was  significantly  lower  in  IDHl-mutant  tumors  compared 
with  those  generated  using  parental  HCT116  cells  (Fig.  4E), 
indicating  that  TCA  metabolism  may  also  be  perturbed  in 
tumors  comprised  of  IDHl-mutant  cells.  In  addition,  the 
contribution  of  glutamine  anaplerosis  to  the  aKG  pool  was 
significantly  elevated  in  IDHl  R132H/-f  and  IDHl  R132C/-t- 
tumors  (Fig.  4F),  which  was  also  observed  in  our  MFA  results 
(Fig.  IH  and  I).  Thus,  the  available  data  are  consistent  with 
our  in  vitro  MFA  results  and  provide  evidence  that  TCA 
metabolism  is  similarly  compromised  by  IDHl  mutations 
in  vivo. 

Inhibition  of  mutant  IDHl  does  not  rescue 
reprogramming  of  TCA  metabolism 

One  possible  explanation  for  the  decrease  in  reductive 
carboxylation  flux  in  IDHl-mutant  cells  is  that  localized 
substrate  (aKG  and  NADPH)  consumption  by  the  mutant 
enzyme  for  production  of  2-HG  compromises  this  activity. 
Therefore,  we  examined  whether  pharmacologic  inhibition 
of  mutant  IDHl  could  increase  reductive  carboxylation 
activity  and  rescue  the  ability  of  ceOs  to  use  this  pathway 
for  growth  under  hypoxia.  To  test  this  hypothesis,  we  treated 
IDHl  R132H/-f  2H1  or  IDHl  R132C/+  2A9  cells  with  a 
mutant  IDHl  inhibitor  (IDHli  A)  similar  to  a  previously 
described  structural  class  (Fig.  5A;  refs.  22,  46).  Doses  of  10 
|lmol/L  were  able  to  decrease  2-HG  levels  more  than  25-fold 
in  both  clones  (Fig.  5B).  As  would  be  expected  from  an 
engineered  cell  line  that  does  not  exhibit  growth  dependence 
on  mutant  IDHl  or  2-HG,  10  pmol/L  of  IDHli  A  had  no 
appreciable  effect  on  the  growth  rate  of  either  cell  line  (Fig. 
5C).  Both  short-term  (3  day)  and  long-term  (31  day)  treat¬ 
ment  with  10  pmol/L  IDHli  A  induced  minimal  changes  in 
metabolite  abundances  beyond  2-HG  (Fig.  5D)  and  effec¬ 
tively  reversed  the  mutant  IDH-dependent  EMT  phenotype 
exhibited  by  these  cells  (Supplementary  Fig.  S5A). 

IDHli  A  did  not  rescue  the  ability  of  cells  to  initiate  reductive 
TCA  metabolism  under  hypoxia,  as  labeling  of  citrate  (Fig.  5E) 
and  other  metabohtes  (Supplementary  Fig.  S5B-S5D)  from 


{Continued.)  F,  aKG/citrate  ratio  from  celis  cultured  as  in  Fig.  5D.  G,  contribution  of  [U-'^Cslglufamine  fo  iipogenic  AcCoA  from  cells  cultured  as  in  Fig.  5D. 

H,  2-HG  levels  in  HCT1 16  IDHl  R132H/+  2H1  and  HCT116  IDHl  R132C/+  2A9  cell  lines  treated  with  the  indicated  concentration  of  IDH1-C227  for  3  days. 

I,  relative  level  of  reductive  glutamine  metabolism,  determined  by  M5  labeling  of  citrate  from  [U-13C5]glutamine  in  HCT116  IDHl  R132H/+  2H1  and 
HCT1 16  IDHl  R132C/+  2A9  cell  lines  treated  with  1  pmol/L  of  IDH1-C227  for  3  or  12  days,  the  final  72  hours  of  which  the  cells  are  grown  in  hypoxia  (2% 
oxygen).  J,  aKG/citrate  ratio  from  cells  cultured  as  in  Fig.  51.  K,  contribution  of  [U-^^Cslglutamine  to  Iipogenic  AcCoA  from  cells  cultured  as  in  Fig.  51. 
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IDH1  Mutations  Sensitize  Cells  to  Mitochondrial  Inhibitors 


[U-^^Cs] glutamine  was  not  increased  compared  with  vehicle 
treatment.  Other  indices  of  reductive  TCA  metabolism,  includ¬ 
ing  the  ratio  of  aKG/citrate  and  contribution  of  glutamine  to 
lipid  biosynthesis,  also  indicated  that  IDHli  A  failed  to  rescue 
reductive  carboxylation  flux  in  these  cells  (Fig.  5F  and  G).  At  10 
|tmol/L,  the  dose  that  showed  maximal  2-HG  inhibition,  IDHli 
A  mildly  inhibited  reductive  carboxylation  in  the  WT  parental 
cells  (Supplementary  Fig.  S5E  and  S5F),  potentially  due  to  off- 
target  effects  on  WT  IDHl  at  high  concentrations.  To  further 
address  this  issue,  we  also  treated  the  IDHl-mutant  cells  with 
an  additional  inhibitor  of  mutant  IDHl  at  more  moderate 
concentrations  (Fig.  5H).  We  again  observed  no  rescue  in 
reductive  glutamine  metabolism  (Fig.  5I-K),  providing  evi¬ 
dence  that  inhibition  of  IDHl-mutant  activity  may  be  insuf¬ 
ficient  to  remove  the  block  in  metabolic  reprogramming  in 
response  to  hypoxic  stress. 

Cells  expressing  mutant  IDHl  are  sensitive  to 
pharmacologic  inhibition  of  mitochondrial  oxidative 
metabolism 

In  comparing  the  growth  rates  of  the  HCT116  panel  under 
normoxia  and  hypoxia,  we  observed  that  mutant  IDHl  cells 
grew  more  poorly  under  conditions  of  low  oxygen  tension 
than  parental  cells  or  those  expressing  mutant  IDH2  (Fig. 
6A).  HCT116  IDHl-mutant  xenografts  also  grew  at  a  signif¬ 
icantly  slower  rate  than  the  HCT116  parental  cells  (Fig.  6B), 
conditions  that  exhibited  significant  stabilization  of  HIFla 
in  both  parental  and  IDHl-mutant  tumors  (Fig.  6C).  The 
growth  rate  of  IDH2-mutant  cells  as  xenografts  was  also 
significantly  decreased  relative  to  the  parental  cells  (Sup¬ 
plementary  Fig.  S6A),  though  our  available  data  do  not  yet 
support  a  hypothesis  for  how  IDH2  mutations  affect  in  vivo 
growth. 

The  slow  growth  of  IDHl-mutant  cells  in  the  xenograft 
model  suggests  that  altered  TCA  metabolism  may  contribute 
to  the  slower  growth  of  these  cells  under  the  decreased 
oxygen  levels  in  vivo.  IDHl-mutant  cells  also  exhibited 
increased  oxidative  TCA  metabolism  under  hypoxia  com¬ 
pared  with  parental  cells  (Fig.  IE  and  I),  providing  evidence 
that  they  are  more  dependent  on  OXPHOS.  To  further 
confirm  this  phenotype,  we  measured  oxygen  consumption 
in  parental  and  IDHl-mutant  cells  under  both  normoxia  and 
hypoxia.  Consistent  with  our  MFA  results,  basal  respiration 
was  not  significantly  altered  in  parental  and  IDHl-mutant 
cells  under  normoxic  conditions,  though  uncoupled  respi¬ 
ration  was  decreased  in  the  IDHl-mutant  cells  (Supplemen¬ 
tary  Fig.  S6B-S6E).  Notably,  mutant  IDHl  cells  exhibited 
significantly  higher  oligomycin-sensitive  oxygen  consump¬ 
tion  under  hypoxia  compared  with  parental  HCT116  cells,  an 


effect  not  reproduced  under  normoxia  (Fig.  6D  and  E  and 
Supplementary  Fig.  S6F  and  S6G).  Therefore,  we  hypothe¬ 
sized  that,  as  with  growth  in  hypoxia,  cells  harboring  IDHl 
mutations  may  be  more  susceptible  to  inhibition  of  oxida¬ 
tive  mitochondrial  metabolism  than  cells  with  WT  IDHl/2 
or  mutant  IDH2. 

To  address  this  question,  we  cultured  parental  HCT116  cells 
and  three  IDHl-mutant  clones  in  the  presence  of  several 
compounds  that  inhibit  Complex  I  of  the  ETC  and  OXPHOS. 
Confluency  measurements  were  taken  every  12  hours,  and  the 
maximum  specific  growth  rate  of  each  cell  was  determined 
using  a  generalized  logistic  growth  model  and  compared  with 
vehicle  treatment  for  each  cell  line  (Fig.  6F).  The  proliferation 
rate  of  cells  with  mutant  IDHl  was  significantly  more  affected 
than  that  of  parental  HCT116  cells  in  response  to  discreet  dose 
ranges  of  Complex  I  inhibitors.  On  the  other  hand,  the  IDH2 
R172K/+  cells  displayed  no  such  increased  sensitivity  with  the 
same  treatments  (Supplementary  Fig.  S6H).  This  altered  sen¬ 
sitivity  is  not  due  to  differences  in  target  modulation,  as  100 
nmol/L  rotenone  effectively  shut  down  oxidative  TCA  cycle 
flux  in  all  cells  tested  (Supplementary  Fig.  S6I).  Treatment  of 
parental  HCT116  cells  with  100  nmol/L  rotenone  also  induced 
reductive  carboxylation,  whereas  R132H/-f  2H1  and  R132C/-f 
2A9  HCT116  cells  were  less  able  to  increase  flux  through  this 
pathway  (Supplementary  Fig.  S6J  and  S6K).  Treatment  with 
Antimycin  A,  an  inhibitor  of  Complex  III  of  the  ETC,  also 
inhibited  oxidative  TCA  metabolism  (Supplementary  Fig.  S6L). 
Induction  of  reductive  carboxylation  by  Antimycin  A  was 
observed  in  the  parental,  but  not  IDHl-mutant,  cells  (Supple¬ 
mentary  Fig.  S6M  and  S6N).  However,  this  compound  had 
differential  effects  on  succinate  labeling  compared  with  rote¬ 
none  (Supplementary  Fig.  S60).  Notably,  IDHl-mutant  cells 
did  not  exhibit  increased  sensitivity  to  Antimycin  A  (Supple¬ 
mentary  Fig.  S6P),  suggesting  that  Complex  III  inhibition 
suppresses  growth  through  distinct  mechanisms  compared 
with  Complex  I  inhibitors  (e.g.,  reactive  oxygen  species  gen¬ 
eration,  and  pyrimidine  synthesis;  ref  43).  Thus,  these  data 
indicate  that  IDHl-,  but  not  IDH2-,  mutant  cells  are  selectively 
sensitive  to  Complex  I  inhibitors. 

To  determine  whether  IDHl  mutants  are  generally  more 
sensitive  to  other  treatments,  we  examined  the  effect  of  the 
cell-cycle  inhibitor  flavopyridol  in  the  HCT116  panel  of  cells. 
IDHl-mutant  cells  did  not  display  increased  sensitivity  in 
comparison  with  parental  or  IDH2-mutant  cells  (Supplemen¬ 
tary  Fig.  S6Q),  further  suggesting  that  the  differential  sensitivity 
we  observe  is  specific  to  inhibitors  of  mitochondrial  metabo¬ 
lism.  These  results  indicate  that  oncogenic  IDHl  mutations 
induce  cells  to  rely  more  heavily  on  Complex  I  of  the  ETC, 
rendering  these  cancer  cells  more  susceptible  to  inhibition  of 


Figure  6.  Cells  expressing  mutant  IDHl  are  sensitive  to  pharmacologic  inhibition  of  mitochondrial  oxidative  metabolism.  A,  doubling  times  of  HCT1 16  cells 
cultured  in  normoxia  or  hypoxia  (2%  oxygen)  for  72  hours.  B,  growth  curves  for  HCT116  parental  and  IDH1-mutant  xenografts.  C,  Western  blot  showing 
HIFla  expression  in  HCT1 16  parental  and  HCT1 16  IDHl  R132H/+  2H1  cells  grown  in  normoxia  in  cell  culture  (last  two  lanes)  or  as  xenografts.  D,  ATP-linked 
oxygen  consumption  for  the  indicated  cell  lines  grown  in  normoxia.  E,  ATP-linked  oxygen  consumption  for  the  indicated  cell  lines  grown  in  hypoxia  (3%  02). 

F,  growth  charts  from  cells  cultured  as  indicated.  Images  were  acquired  every  12  hours  to  measure  confluency.  Change  in  growth  relative  to  DMSO 
treatment  (AX%)  was  calculated  using  a  generalized  logistics  growth  model  for  batch  culture  and  represents  the  change  in  the  specific  growth  rate  relative  to 
the  DMSO  treatment  for  the  indicated  cell  line.  G,  heatmap  displaying  IC50  values  to  four  inhibitors  of  mitochondrial  metabolism  for  more  than  500  cancer  cell 
lines;  HT-1 080  and  SW1 353  cells  are  indicated.  H,  growth  of  HT-1 080  and  SW1 353  cells  under  the  indicated  concentration  of  phenformin  for  48  hours. 
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Figure  7.  Mutant  IDH1  sensitizes  celis  to  inhibition  of  oxidative  mitochondriai  metaboiism.  Left,  under  normai  growth  conditions,  glucose  is  metaboiized 
oxidatively  in  the  mitochondria,  and  AcCoA  and  iipids  are  derived  mainiy  from  giucose  carbons.  Middie,  in  IDH1  WT  celis,  inhibition  of  oxidative  mitochondrial 
metabolism  (induced  by  growth  in  hypoxia  or  pharmacologic  inhibitors  of  the  ETC)  limits  glucose  flux  to  the  mitochondria,  and  cells  instead  rely  on 
reductive  glutamine  metabolism  via  IDH1  to  provide  carbons  for  AcCoA  generation  and  lipid  synthesis.  Right,  when  oxidative  mitochondrial  metabolism  is 
inhibited,  cells  with  a  mutant  IDH1  allele  are  unable  to  fully  induce  reductive  glutamine  metabolism  and  are  thus  compromised  for  AcCoA  and  lipid  production, 
leading  to  decreased  cell  growth. 


this  pathway  compared  with  cells  with  WT  IDHl/2  or  mutant 
IDH2  alleles. 

Finally,  to  evaluate  whether  these  results  are  relevant  to  other 
cells,  we  interrogated  the  Cancer  Cell  Line  Encyclopedia  (CCLE; 
ref.  32),  which  contains  compoimd  sensitivity  data  across  more 
than  500  cell  lines  for  four  ETC  inhibitors  (Fig.  6G;  ref  47).  Cell 
lines  clustered  well  into  sensitive  and  insensitive  groups,  sug¬ 
gesting  these  compounds  show  consistent  behavior  across  a 
wide  panel  of  cell  lines.  Notably,  IDHl-mutant  HT1080  cells  fell 
into  the  sensitive  group,  whereas  IDH2-mutant  SW1353  cells  fell 
into  the  insensitive  group  (Fig.  6G).  The  differential  effects  of 
mitochondrial  metabolism  inhibitors  were  not  likely  due  to 
HT1080  being  generally  more  sensitive  to  compound  treat¬ 
ments,  as  the  sensitivities  of  HT1080  and  SW1353  to  a  broad 
array  of  more  than  1,300  compoimds  were  within  one  SD  of  each 
other  (Supplementary  Fig.  S6R).  Furthermore,  HT1080  cells 
were  significantly  more  sensitive  than  SW1353  cells  to  phen- 
formin  treatment,  a  compound  that  was  not  included  in  the 
CCLE  screening  set  (Fig.  6H).  Together,  these  data  indicate  that 
IDHl  mutation  may  substantially  sensitize  cells  to  inducers  of 
mitochondrial  stress. 

Discussion 

Since  the  discovery  of  oncogenic  mutations  in  IDHl  and 
IDH2,  significant  efforts  have  been  made  to  elucidate  the 
mechanisms  driving  tumorigenesis  in  these  cancers.  Owing 
to  the  accumulation  of  D-2-HG  in  these  tumors,  researchers 
have  focused  on  the  role  of  this  oncometabolite  in  regulating 
the  phenotype  of  IDHl/2-mutant  cancer  cells.  For  example, 
high  D-2-HG  levels  and  other  metabolites  regulate  the  activity 
of  aKG-dependent  dioxygenases  that  control  many  distinct 
cellular  processes  (12).  However,  the  diverse  roles  of  these 
enzymes  in  mediating  activities  ranging  from  collagen  hydrox- 
ylation  and  HIF  stabilization  to  epigenetics  regulation  com¬ 


plicate  identification  of  the  specific  process(es)  driving  tumor¬ 
igenesis  in  each  tumor  type. 

Despite  the  central  role  of  these  enzymes  in  cellular  metab¬ 
olism,  surprisingly  few  investigations  have  addressed  the  met¬ 
abolic  changes  that  occur  as  a  result  of  these  genetic  mod¬ 
ifications.  Here,  we  find  that  IDHl  mutations  cause  cells  to 
increase  flux  through  the  oxidative  TCA  cycle,  increase  respi¬ 
ration,  and  compromise  the  conversion  of  glutamine  to  citrate, 
AcCoA,  and  fatty  acids  under  hypoxia  compared  with  those 
with  WT  IDHl  (Fig.  7).  Others  have  previously  shown  that 
IDHl-mutant  proteins  are  biochemically  compromised  with 
respect  to  this  latter  functionality,  suggesting  that  cells  har¬ 
boring  such  mutations  may  be  similarly  defective  under  certain 
conditions  (26).  However,  the  cellular  consequences  of  this 
effect  have  not  been  well  characterized  within  intact,  hetero¬ 
zygous,  IDH-mutant  cells.  Proliferating  cells  must  double  their 
membrane  lipids  to  successfully  complete  cell  division,  and 
evidence  suggests  that  tumors  may  rely  more  on  de  novo 
lipogenesis  than  do  nonneoplastic  tissues,  and  inhibition  of 
lipid  synthesis  decreases  tumor  growth  in  vivo  (48,  49).  In 
addition,  AcCoA  is  an  important  precursor  for  a  number  of 
other  molecules,  including  cholesterol,  phospholipids,  amino 
acid  modifications,  and  histone  acetylation  (50).  Interestingly, 
previous  studies  have  also  found  that  overexpression  of  mutant 
IDH  leads  to  a  decrease  in  N-acetyl  amino  acids,  and  these 
changes  were  also  observed  when  comparing  WT  human 
glioma  tissue  with  that  of  tumors  with  mutant  IDHl  (24).  This 
suggests  that  other  AcCoA-dependent  molecules  may  be  sim¬ 
ilarly  perturbed  in  the  IDHl-mutant  setting.  Our  application  of 
MFA  to  IDHl-mutant  cells  builds  upon  these  results  by  addres¬ 
sing  the  functional  consequences  of  heterozygous  IDHl  muta¬ 
tions  and  in  particular  the  metabolic  limitations  that  arise  in 
tumors  cells  with  these  genetic  modifications.  Given  the 
importance  of  each  of  these  AcCoA-dependent  processes  for 
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cellular  homeostasis  and  proliferation,  we  speculate  that  the 
reduced  metabolic  flexibility  of  these  cells  contributes  to  the 
decrease  in  growth  rate  that  we  observe  in  the  IDHl-mutant 
cells  under  conditions  of  decreased  oxidative  mitochondrial 
metabolism. 

Although  2-HG-mediated  control  of  aKG-dependent  diox¬ 
ygenase  activity  clearly  plays  a  role  in  tumorigenesis  driven  by 
IDH  mutations  (13-18, 20-23),  our  results  provide  insights  into 
therapeutic  strategies  that  exploit  the  metabolic  vulnerabilities 
caused  by  partial  loss  of  WT  IDHl  function.  Interestingly,  we 
observe  that  IDHl-mutant  cells  do  not  exhibit  pronounced 
metabolic  differences  in  normoxia;  however,  growth  in  low 
oxygen  tensions  or  with  pharmacologic  inhibitors  of  mito¬ 
chondrial  metabolism  results  in  the  emergence  of  dramatic 
metabolic  changes.  Our  MFA  results  identify  several  enzymes 
and  pathways  that  are  altered  under  hypoxia  and  in  particular 
in  the  context  of  IDHl  mutations.  Although  compartment- 
specific  IDH  fluxes  cannot  be  resolved  with  these  data,  these 
findings  further  highlight  the  importance  of  WT  IDHl  activity 
in  mediating  reductive  glutamine  metabolism.  Importantly, 
our  results  provide  evidence  that  IDHl  mutations  functionally 
compromise  cellular  metabolism  under  conditions  of  low 
oxygen  levels,  with  the  most  pronounced  effects  being 
increased  dependence  on  oxidative  mitochondrial  metabolism 
and  an  inability  to  induce  reductive  glutamine  metabolism.  We 
artificially  induced  such  stresses  using  pharmacologic  inhibi¬ 
tors  of  Complex  I  or  manipulation  of  the  oxygen  tensions,  and 
observed  selective  growth  rate  reductions  in  several  IDHl- 
mutant  cells,  but  not  in  parental  or  IDH2-mutant  cancer  cells. 
Other  recent  studies  have  also  highlighted  the  importance  of 
oxidative  mitochondrial  metabolism  for  tumor  cell  growth  and 
survival  both  in  vitro  and  in  vivo  (45,  51). 

These  results  suggest  that  compromised  IDHl  function  may 
affect  the  proliferative  capacity  of  tumor  cells  and  furthermore 
that  IDHl-mutant  tumors  may  be  sensitive  to  inhibitors  that 
perturb  mitochondrial  metabolism.  When  comparing  the  met¬ 
abolic  phenotype  of  tumor  xenografts  derived  from  parental  or 
IDHl-mutant  cells  with  our  in  vitro  results,  similar  changes 
were  detected,  including  increased  glutamine  anaplerosis  and 
a  decreased  aKG  to  citrate  ratio.  The  increase  in  glutamine 
anaplerosis  we  observe  in  the  IDHl-mutant  cells  is  in  agree¬ 
ment  with  previous  findings,  which  suggest  that  IDHl-mutant 
cells  display  an  increased  sensitivity  to  glutaminase  inhibitors 
(25).  The  similar  metabolic  changes  that  could  be  rehably 
measured  in  vivo  suggest  that  the  altered  sensitivity  we  observe 
to  inhibitors  of  mitochondrial  metabolism  in  vitro  may  also  be 
true  in  vivo.  Additional  studies  are  required  to  determine  if 
cellular  proliferation  in  the  tumor  microenvironment  alone 
can  drive  hypoxia  and  induce  reductive  glutamine  metabolism. 
Regardless,  tumors  would  stiU  be  expected  to  increase  their 
reliance  on  WT  IDHl  (or  cytosolic  TCA)  activity  when  treated 
with  phenformin  or  other  inhibitors  of  mitochondrial  metab¬ 
olism,  suggesting  that  these  strategies  could  be  efficacious  in 
IDHl-mutant  cancers.  As  such,  this  increased  susceptibility  of 
cultured  IDHl-mutant  cells  relative  to  parental  cells  or  IDH2- 
mutant  cells  provides  intriguing  evidence  of  a  potential  ther¬ 
apeutic  strategy  associated  with  IDHl  mutational  status  and 
warrants  further  investigation  in  preclinical  models. 


We  find  that  inhibition  of  mutant  IDHl  is  unable  to  reverse 
the  observed  metabolic  phenotype.  The  DNA  hypermethylator 
phenotype,  which  is  highly  associated  with  IDH  mutation,  is 
also  not  entirely  reverted  by  a  mutant  IDHl  inhibitor  (22), 
providing  further  evidence  that  some,  but  not  all,  mutant  IDH- 
dependent  phenotypes  may  be  reversed  by  inhibitors  targeting 
2-HG  production.  Mechanistically,  this  result  also  suggests  that 
the  metabolic  defect  we  observe  may  be  independent  of  2-HG 
production.  A  previous  study  used  biochemical  assays  to 
quantify  the  effects  of  IDHl  mutations  on  reductive  carbox- 
ylation  activity,  and,  in  agreement  with  our  findings  here, 
demonstrated  that  the  mutant  enzymes  are  unable  to  catalyze 
the  conversion  of  aKG  and  CO2  to  isocitrate  (26).  This  study 
concluded  that  the  subunits  in  a  WT/mutant  heterodimer 
function  independently;  however,  our  modeling  data  indicate 
that  heterozygous  IDHl  mutations  lead  to  a  much  greater  than 
50%  inhibition  of  reductive  glutamine  metabolism  (Fig.  IF), 
suggesting  a  possible  dominant  effect  of  the  mutant  protein  in 
cells  or  alternatively  global  metabolic  reprogramming  in 
response  to  the  compromised  cytosolic  IDHl  activity  caused 
by  these  mutations. 

Importantly,  as  small  molecules  capable  of  inhibiting 
mutant  IDHl  enzymatic  activity  and  preventing  D-2-HG 
accumulation  fail  to  rescue  mutant  cell  metabolism  under 
hypoxia,  this  suggests  that  combinatorial  therapeutic  strat¬ 
egies  that  block  oncogenic  D-2-HG  production  (e.g.,  via  a 
mutant-selective  inhibitor  of  enzyme  function),  while  simul¬ 
taneously  targeting  mutant  IDHl-induced  metabolic  liabil¬ 
ities,  may  be  a  viable  option  for  therapy.  Such  an  approach 
could  involve  IDHl-mutant  inhibitor  treatment  to  attenuate 
any  prosurvival  or  dedifferentiation  effects  of  D-2-HG,  while 
increasing  the  tumor's  reliance  on  WT  IDHl  activity  through 
an  inhibitor  of  oxidative  mitochondrial  metabolism.  As 
drugs  that  could  target  the  mutant  IDHl  metabolic  pheno¬ 
type  are  already  in  the  clinic  (metformin,  phenformin)  and 
inhibitors  of  mutant  IDHl  are  currently  being  developed 
(22,  46),  it  is  hoped  that  this  hypothesis  will  be  tested  in  the 
clinic  in  the  near  future. 
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